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ABSTRACT 


A model of the elongation characteristics of a double 
braid nylon rope is presented as an aid to sizing a marine 
towline. The material and structural elongation properties 
are combined and modeled as being the sum of a non-elastic 
permanent elongation and a load dependent elastic elongation. 
The loading range of practical interest is examined and a 
design factor of safety, with both an upper and lower bound, 
is established. The lower bound 1s incorporated to limit 
relative movement between structural elements in an attempt 
to control internal abrasion. The upper bound is imposed to 
stay within the safe working load and total elongation limits 
of the rope. The governing equations for the towline are 
then solved within these constraints for a typical submarine 
towing system to provide a towline capable of sustaining the 
hydrodynamic loads due to the resistances of the submarine 
and the towline. It is shown that the resistance of a 
typical towline may be of the same order of magnitude as the 
resistance of the towed vessel and cannot be neglected. The 
solution results in a recommended towline diameter and break- 
ing strength which is then examined under possible loading 
conditions other than those selected for the initial design. 
It 1s shown that a significant portion of the operating 
envelope is outside the limits obtained from the factor of 
safety analysis and that a nylon towline may be vulnerable 
to significant internal abrasion, especially if the line is 
Oversized. Alternate materials are discussed but not anal- 
yzed due to the lack of experimental data which could sup- 
port a functional model of elastic elongation behavior. 
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enapter I 
INTRODUCTION 


A current problem in ocean engineering is the incorpora- 
rom Of synthetic rope into lifting, handling, mooring, and 
towing systems. The high strength-to-weight ratios, the 
energy absorption characteristics, and the flexibility and 
ease in handling are all features offered by synthetic rope 
that are not common to wire rope or natural fiber rope. The 
adaptation of these characteristics to marine engineering prob- 
lems can result in improved system safety and performance if 
the design is well executed. The design of an efficient ship- 
to-ship towing system requires an analysis which brings 
together those material and structural characteristics needed 
to sustain the static and dynamic loads imposed by the con- 
nected vessels as they respond to the excitation of the ocean 
environment. 

One of the most important design elements of the towing 
system is the ability to absorb dynamic loads and dissipate 
energy without damaging system components or degrading system 
performance. This can be accomplished through employment of 
mechanical tensioning devices, steel/chain catenaries, syn- 
thetic towlines and, in some special cases, the use of the 
control surfaces on the towed vessel. The application of a 
synthetic towline is the most universal of these methods as 
it is easily adapted to both scheduled and emergency tows. 
The resort to synthetic rope on such occasions has been com- 
mon for several years, but has recently been discouraged, and 
even forbidden by the U. S. Navy in scheduled military tows, 
and by Lloyds of London for insured commercial tows. Both of 
these actions were apparently the result of the unacceptably 
high failure rate of synthetic towlines and the resulting 
damage. There are many possible causes of towline failures. 
While it is not the intent of this thesis to address them all, 
we submit that failure rates can be reduced and conditions 
leading to failure can be more accurately predicted by im- 


proved design processes. 
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This thesis will address the incorporation of synthetic 
rope into the marine towing system. Nylon, polyester, poly- 
propylene, and Kevlar are the four primary fibers used to 
make synthetic rope, each with distinct mechanical charac- 
teristics capable of improving the system being designed. 
However, such variations in characteristics also require 
that slightly different methods be used in the analysis. 
Seeteewitl be used here as the base material for the presen- 
tation because it encompasses the vast majority of the ropes 
presently in use. Even though the methods to be suggested 
will be adaptable to other materials, it will be imperative 
that differences in mechanical characteristics be accounted 
for. In addition, the structure of the rope will be limited 
to 2-in-l double braid since it is a torque-free structure 
and offers the highest breaking strength for a given material 
and size when compared to the other common commercial rope 
Seructures. 

This design process will present a method of predicting 
the elastic behavior of a nylon towline, examine the effects 
of the design "factor of safety", present a quasi-static pre- 
liminary design and examine the sensitivity of the design to 
each of the design parameters. This process of examining the 
interaction between material properties, loading conditions, 
and design methodology will lead to the preliminary sizing of 
the synthetic rope and to subsequent improvements of the tow 


system through trade-offs between design parameters. 
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Chapces ie 


ELONGATION CHARACTERISTICS 


mee Background 


To describe the elastic elongation characteristics of 
synthetic ropes requires the selection of the best working 
approximation to a complex mechanical structure comprised 
of a very large number of basic filament elements. Polymeric 
fiber elements are, in fact, viscoelastic in nature; any 
approximations that ignore time dependence are inherently 
inexact. Polymeric fibers subjected to a sustained load, 
such as the mean towing load, are known to creep. This 
creep is reflected in an increase in nonrecoverable elonga- 
tion and the fiber will eventually rupture; this creep rup- 
ture is a function of load amplitude and time under load. 

A fatigue failure model has been developed for predicting 
failure of low twist yarns based on creep rupture data from 
fibers [11], but this model has not yet been extended to 
marine rope level. The inclusion of the time factor in a 
truly viscoelastic model which could be matched with the 
complex time dependent excitation forces of the ocean envi- 
ronment cannot be done at this time. 

The task is further complicated by the need to design 
towing systems requiring the use of large ropes which have 
not been extensively tested. Comparison of data which have 
been collected by various sources must be undertaken with 
caution [4], since test methods, environmental conditions, 
and specimen size and geometry greatly influence results 
obtained. Several ongoing development programs [8] are 
underway, focusing on various aspects of synthetic fiber 
rope behavior in the marine environment. When such programs 
have been completed, it should be possible to eliminate many 
of the unknowns that now hinder the design process. [In the 
interim some uncertainties must be accommodated in the de- 


Sign methods used. One such method is presented here. 
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weep? Load—-Elongation Characteristics 


It is essential that the elongation characteristics 
peculiar to synthetic ropes be understood at the outset. An 
introduction to the subject is here provided with the aid of 
Fig. 1. A more complete description can be found in [1, 2, 
and 3]. 

As a synthetic rope is loaded for the first time from 
T, to a working load T,, the load-elongation curve follows 
the path A-B of Fig. l(a). If the load is then decreased to 
Sepmerve second portion of the curve B-C is traced; as the 
load is again increased to To the path C-D is formed and the 
typical hysteresis loop can be seen. The distance A-C 
represents the non-elastic elongation, AL. part of which 
would be recovered if the rope were allowed to relax for ex- 
tended periods between cycles. Subsequent cycles will cause 
the hysteresis loop to migrate to the right as the non- 
elastic elongation is increased with each load cycle. The 
hysteresis loop "stabilizes" [1] at approximately fifty 
cycles and it is this loop that is used for design applica- 
Miens that require cyclic loading. 

Figure 1(b) shows the first, tenth, and fiftieth cycles 
Siecameyoical rope stabilization profile [1] with the perti- 
nent components of the total elongation labeled. It can be 
seen that the total elongation is composed of the non-elastic 
elongation and the working elongation. These elongations can 
be converted to strains by dividing by a "reference" length 
and at this point extreme care should be taken when comparing 
data from different sources because the same "reference" 
length is not always used. For example, reference [1] de- 
fines working strain as the change in working length over 
the working length, where the working length is the base 
length L, plus the non-elastic elongation AL. In contrast, 
reference [2] defines the working strain as the elastic 


change in length over the base length. The relation is: 
SS L/L (1) 


For nylon working between twenty and thirty percent of the 
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average breaking strength, Bo, L/L, ifseappLeoOximareLy 1.13 
[1]. Another contrast is found in Reference [4], for the 
testing of wet single-point-mooring hawsers, where the ref- 
Sence length was measured at T, with the rope wet. In this 
case, a direct conversion is not possible for nylon because 

of the shrinkage which causes the fiber to swell and the length 
to shorten. However, the author [4] states eiehe telat results 

of a few tests were re-analyzed and that the elasticity re- 


ported was about 1% above that based on the working length. 


Pee Elastic-Elongation Model 


With these definitions it is proposed [1] that rope size 
can be eliminated from the stress-strain relation by dividing 
the actual tension, T, by the average rated breaking strength, 
Ba, to produce the specific tension, tT, which can then be re- 


lated to the working strain by a function of the form: 
= m 
t = Al x (e,) (2) 
where Al and m are empirically determined constants and the 
rope has been stabilized with 50 load cycles at 20% of the 


breaking strength. For nylon in both the wet and dry condi- 


tion, the function becomes: 


aT 9.78 x (emi = (wet nylon) (3) 


a 
Il 


Arex (ile (dry nylon) (4) 


In addition to these functions, Reference [5] presents 
a design curve for a 2.5-inch double braid rope subjected to 
cyclic load conditions. Again, the design curve is based on 
eeeaecOllected during the fifty-first cycle at 20% of Be. 
Reference [6] presents a "linearized (retention) equation" 
for nylon towlines which is size dependent; this function is 
presented as part of a towing system design report and the 
basis for the equation is not given. This report uses E for 


the elastic elongation of the towline in feet, and is: 
Ee (ceca my (7.4 xO" x d*) (5) 


Since E is stated to be the elastic elongation, it is 


Pre og 





assumed to be equivalent to AL. By rearranging the terms and 
dividing the tension with the rated breaking strength, to 
produce a function of specific tension for the given size rope, 
the equation can be compared with the above functions in the 


Zmorm: 


Baby = (T/B_)/1(7.4 x 10” & a 31 (6) 


where it has been assumed that the length intended was the 
meeeeengthn, L,. The results of this equation can now be 
plotted for ropes having diameters in the range of interest 
and compared to the size independent functions, equations (3) 
and (4). 

Mier. 92, a Common “reference” length of L, has been 
established and for the functions given by McKenna [1], the 
data given by Flory [4], the typical design curve presented 
by Wong [5], and the linear function used in the towing 
system report [6], the respective elastic elongation curves 
are presented for approximate comparison. In Fig. 2a, curve 
(1) represents the wet nylon function, WNF, and curve (2) 
represents the dry nylon function, DNF, both of which are size 
independent, but based on data taken after the fiftieth load 
cycle to 202 Be: Curve (3) represents the design curve, DNC, 
given by Wong for 2.5-inch diameter dry double braid nylon 
meeenwon the fifty-first cycle. Curves (L1.75) and (ZL 2.5) 
represent the "linearized (retention) equation" for 1.75- and 
2.5-inch diameter ropes. It appears that the linear function 
is fairly accurate at low specific loads when compared to the 
evemy on function and the dry nylon curve. [In Fig. 2b, 
curves (4) through (9) represent data taken at 10, 100, 300, 
1,000, 3,000, and 10,000 cycles [4] for 1.75-inch diameter 
pre-soaked double braid nylon ropes. These curves for pre- 
soaked double braid nylon ropes more closely follow the dry 
nylon function. The reaction of the pre-soaked ropes can be 
partially justified by assuming that the water is being 
Squeezed out as the rope is cycled under load which results 
in increased friction between structural elements. This 


would probably produce only small changes in the elastic 
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behavior of the rope. It can also be assumed that the rope 
was beginning to dry out [4] from internal heat generation and 
exposure to the air during testing; approximately 2,500 cycles 
were applied each day [4]. Both of these factors would tend 
to shift the curves toward the dry nylon model, as they are 
shown in Figs. 2a and 3. In Fig. 3, curve (3) was eliminated 
because of the uncertainty in the original reference length 
and only the curves (5) and (8), representing 100 and 3,000 
cycles, were retained from the pre-soaked cyclic load tests. 
Once again these are compared to the linear functions 
presently used to design some U. S. Navy tow systems. 

Figure 4 presents the non-dimensional stiffness, i.e. 
the change in percent load with a change in percent elonga- 
tion, and it is here that the validity of the linear approxi- 
mations becomes gquestionable. If it is assumed that the 
curve (1) for wet nylon is closest to the actual rope behav- 
ior, then the linear approximation is valid only in the 
region of the 15% load range and that there is a deviation of 
between 30% and 50% when the linear function is compared to 
Bieewet nylon function, curve (1), for a cyclic load ampli- 
tude about the mean load of plus or minus 10%, as is evident 
when the respective values are compared at tT equal to 5% and 
25%. Such a load amplitude could reasonably be expected under 
working conditions. In contrast, if it is assumed that the 
meemyton function, curve (2), or the 100 and 3,000 cycle 
curves (5) and (8) approximate the real condition, then the 
linear approximation is accurate only very near a specific 
load of 5%, which should not represent normal working con- 
ditions. 

Based on these comparisons no further consideration will 
be given to the linear approximation of the load-elongation 
behavior of nylon towlines, since it is apparent that both 
the elongation and the stiffness are very non-linear below 
Bos of B, and it is in this region that we are most interested. 

To this point, it has been assumed that a rope is com- 
pletely stabilized by the fiftieth load cycle and that 
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further creep migration of the hysteresis loop will not be 
Significant, except in the case where the rope is unloaded 
and allowed to recover slowly for an extended period, ch 
which case the rope may need to be "broken-in" again. The 
fiftieth cycle stability assumption is the basis for Eqs. 

(3) and (4) and for the proposed design curve [5], curve (3). 
In fact, the assumption is based on synthetic rope test 
Standards [7]. This is not a true representation for an 
extremely high number of load cycles that could be expected 
during a towing operation (5,000 to 15,000 cycles per day) 
which may encompass several days of continuous operations 
where there iS a variable but constantly applied mean load. 
The hysteresis loop actually has a nearly logarithmic mi- 
gration [1], with the exceptions noted above, and the perma- 
nent elongation will continue to increase as the elastic 
elongation decreases. This iS simply the manifestation of 
the viscoelastic properties of the nylon filaments and even 
though it cannot be quantified in the complex structure of a 
large rope, it must be accounted for. Referring to Fig. 2a, 
it can be seen that such a migration would cause the high 
cycle wet nylon curve (1) to shift again toward the low 
cycle dry nylon curve (2). Similar to the shift seen in 
Fig. 2b, as the number of cycles is increased from 10 cycles, 


curve (4) to 10,000 cycles curve (9). 


2.4 Comparison with Data 
For these reasons, the elastic elongation function for 
dry nylon, Eq. (4), will be used in this thesis, recognizing 
eiat it is: 
(1) Approximate, 
(2) Chosen to represent high cycle wet behavior. 
Under conditions where a new towline was being installed 
for a specific task, it might be prudent to use the wet 
nylon function to predict initial load elongation behavior, 
providing a lower bound on load amplification and an upper 
bound on elongation. Then, to use the dry nylon curve to 


provide the upper bound on load amplification and lower bound 
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on elongation, the condition that could be expected after the 
towline had undergone a significant number of cycles. 

The basic assumption thus far has been that the load 
elongation behavior can be normalized by dividing by the 
average rated breaking strength to eliminate the effects of 
changes in the size of the rope. This assumption was checked 
by comparing the predicted stiffness, using the function shown 
in Fig. 4, with an apparent spring constant, Kap, presented by 
Bitting [3]. This was done by dividing the apparent spring 
constant by the manufacturers published breaking strength to 
produce a specific apparent spring constant which is a func- 


tion of frequency, £, mean load, Ta! and load amplitude, DT: 
Kaps = Kap/B. (We 


Bitting presents coefficients to be used in a Box- 
Behnken equation (7a) for the apparent spring constant for 
faeoy 4, 1", and 1%" double braid nylon line: 


Pte ee rose ime B44 x pT= + 


Kap 
+ B6 x Tm + B7 x DT + BS x £ x Tm + BI x 


Shoe aed ex Timex DT C7 a tee 


All of these coefficients were obtained by testing the ropes 
in the wet condition and are considered to be valid only in 
Pmemrange Of testing [3], which is for a specific load greater 
than 10%. The comparison is shown in Fig. 5, where it is 
seen that in the 20% load range the dry nylon function pro- 
Vides a good model for Bitting's %", 1", and 1)" data. The 
departure of the 3/4" apparent spring constant cannot be ex- 
plained at this time. Based on the relative magnitude of the 
coefficients in Eq. (7a), it is found that the apparent 
meemig Constant, Kap, iS primarily a function of mean load, 
Tm, and load amplitude, DT. The frequency of loading, f, 

has a minimal effect within the range of interest. In Fig. 
2a the maximum permissible frequency and load amplitude as 


given by Bitting were used to plot the curves, while in Fig. 
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5b, the minimum permissible values, also given by Bitting, 
memgemused. A comparison of these two figures reflects the 
minimal effect of varying the frequency and load amplitude 
over the total valid range. 

Further load elongation data were presented by Gibson 
and Wolfe [12] for 4%" circumference double braid nylon rope 
tested in both the dry and wet condition. These data are 
presented in Figs. 6 and 7, where the wet and dry nylon func- 
tions have been added for comparison. It is interesting to 
note that these curves represent tensile tests where the 
specimens had been conditioned with 250 cycles to the base 
Bemerence load, T,, and then loaded in a tensile break test. 
This base reference load 1s a very small percentage of the 
breaking strength, but the elastic elongation in the lower 
load ranges remains in good agreement with the proposed 
model. From Fig. 6, it can be seen that the dry tensile 
specimen falls between the dry and wet functional models; 
further, the test results are very close to the dry-nylon 
functional model throughout the range of practical working 
loads. The wet tensile specimen shown in Fig. 7 is very 
close to the wet-nylon functional model in the lower load 
range and indicates good agreement with the lower bound 
stiffness referred to earlier. A second point that must be 
mentioned here is that the ropes tested by Gibson and Wolfe 
were manufactured by Wall Rope Works; this is significant in 
that all of the previous data were taken on ropes manufac- 
tured by Samson, and it was believed that there could be 
variations between different manufacturing companies. This 
concern may still be valid, but based on this admittedly 
limited comparison, the variations may not be a major factor 
in the preliminary sizing phase of the cCesign. 

It should be noted at this point that polyester and 
polypropylene reaction with water at the molecular level is 
negligible, at temperatures commonly found in the ocean en- 
vironment, and thus do not exhibit significant shrinkage or 
strength reduction when exposed to the marine environment. 


A single elastic elongation function may be sufficient to 
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describe their behavior [1]. These functions are based on 
the stabilized rope and caution should be used when applying 
them directly to very high cycle loading conditions because 
of the possible continued creep migration of the hysteresis 
loop. Test data on large polyester and polypropylene is not 


weeeavallable to attempt to verify either function. 


20 = 


Gimbeer fit 


FACTOR OF SAFETY 


me)6DeCfinition 


Modeling elastic elongation as a function of specific 
loading which requires that the working load be approximated 
and that the breaking strength of the rope be known, implies 
that the size of the rope is known apriori. This is not the 
case, since the objective of the design is to size properly 
the rope used in making the towline. An alternate approach 
is to obtain a nominal specific loading by defining an accep- 
table factor of safety. Defining the ratio of the breaking 
strength to the working load as the factor of safety, FS: 


FS = B./T,, (8) 


which for nylon must be further modified to account for the 
reduction in breaking strength due to immersion in the marine 
environment. A range of 10% [3,4] to 20% [8,4] of breaking 
strength 1S commonly accepted with a value of 15% commonly 
melt, 2, 6]. Using 15% here, the Marine Factor of Safety 


for nylon is defined as: 
MFSN = 0.85 x FS (9) 
This implies that the specific loading will be: 


t = 1/FS = 0.85/MFSN (HO). 


3-2 Range 


Factor of safety values range from 1.5 [9], for single- 
point-mooring hawsers, to 9 [10] for general applications. 
For double braid nylon rope being used in high-cycle towing 
applications, there are several considerations which can be 
used to narrow this range. 

It has been found that tensile fatigue failure (due to 
creep) is not so great a problem when ropes are cycled below 


30% of breaking strength, if the minimum load is allowed to 
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go to zero during each cycle. It has also been shown that 
the rope can be cycled as high as 60% of breaking strength 

if the minimum load is never allowed to reach zero and will 
survive a comparable number of cycles. A load of 3% of 
breaking strength is cited [3,9] as an acceptable lower limit 
in this case. McKenna [{1] recommends that the maximum load 
never exceed 30% of breaking strength for safety and that 
cyclic loads not exceed 20% of breaking strength. 

There are also indications that under fast dynamic load- 
ing conditions, the rope will stiffen and have an apparent 
modulus several times that obtained from moderate cyclic load 
tests, magnification factors of 3 to 4 are cited by Bitting 
[3]. While McKenna [1] states, "Loading rates should not ex- 
ceed 5% of the breaking strength per second with a delay of 
no less than ten seconds between cycles to approach a normal 
response", in towing applications, the delay between cycles 
will be dictated by sea state and heading, and it may not be 
possible to meet the ten-second recommended minimum, while a 
properly sized towline could remain below a load rate of 5% 
of breaking strength per second augmented by adjusting the 
speed of the towing vessel. In the event of the occurrence 
of a fast-dynamic load, the larger diameter rope, having a 
higher breaking strength and greater stiffness, would signi- 
ficantly increase the load transmitted to the other compo- 
nents of the system if a dynamic stiffness of 3 to 4 times 
were realized. Stiffness is also known to increase with 
extended exposure to water; this increase can be as high as 
99% [3] with a concurrent strength loss of 50%. 

The towing hawser will be a substantial investment and 
those used by the U. S. Coast Guard are expected to last for 
up to five years, during which time it is often exposed to 
water. The very nature of marine towing is one of cyclic 
loading and the possibility of minimum loads that approach 
zero cannot be excluded (such as a low speed tow in a high 
sea state). The combination of these constraints imply that 
loads should be kept below 30% of breaking strength, but that 


the minimum design load should be greater than 3% of breaking 


1G 


strength range to limit the possibility of zero minimum 
cyclic load conditions. Thus, by examining the behavior of 
double braid nylon rope, a range for the factor of safety and 


the acceptable specific loading has been established: 


<= FS <= 30 (11) 
3.5 <= MFSN <= 35 (12) 
3% <= T <= 29% C13)) 


Note that the lower limit 1s maintained to account for 
the apparent effects of internal abrasion, caused by relative 
motion between the structural elements of the rope. Main- 
taining a minimum load has the effect of tightening the total 
structure and reducing the structural relaxation which results 


in the major component of relative motion. 


meeeeDynamic Load Amplification Factor 


The upper limit, as stated, must include both the static 
and dynamic loads. In a quasi-static preliminary design the 
dynamic loads can be accounted for through the use of an 
assumed dynamic amplification factor. The magnitude of this 
dynamic amplification factor must be selected for each indi- 
vidual towing configuration based on an analysis of the 
relative accelerations between the towing and the towed 
vessels. For a towing system where the towed body 1s a sub- 
merged submarine, which is not excited by waves, the only 
source of dynamic loading is from the towing vessel. Under 
these conditions Kline and Blockwick [6] use a dynamic am- 
plification factor of two which is based on towship motions 
obtained from scale model testing. 

The range for T can now be used in conjunction with the 
elastic-elongation function, curve (2) of Fig. 3, with pre- 
dicted loads and displacements to size properly the required 
towline in diameter and length. This is done by modifying 
the upper limit by dividing by the appropriate dynamic am- 
plification factor, which will be taken as two for the pur- 


pose of this presentation: 


3% <= Tt <= 15% is al) s. 


a= 








Chapter Pays 
QUASI-STATIC MODEL 


4.1 Towing Configuration 


A marine towing evolution in the presence of waves is 
obviously a dynamic environment and an exact evaluation of 
the loads imposed on the components of the towing system 
would include a combination of the static loads, due to the 
steady motion of the towed vessel through the water, and the 
dynamic loads of the time-varying wave excitation. Dynamic 
analysis of the towing system 1S an area of current research 
and is beyond the scope of this presentation. A quasi-static 
approach will be used here, with the application of the dy- 
namic amplification factor imposed on the specific loading to 
account for wave excitation. 

The model will be developed around an assumed towline 
profile for a submerged towing operation similar to that shown 
in Fig. 8. The towed vessel, at Point A, will be modeled as 


a notional submarine having dimensions as shown in Table l. 


TABLE 1 


Design Point Geometry and Constraints 


Dimensions Design Constraints 


Length 200 ft. Maximum Submerged Towing Ve- 
iYeyen eng 15 knots 


Beam 25 ft. Design Tow Depth Z0CR ie 
Maximum Vertical Tension 
ae 15,000 Bes. 
av 
Nominal Vertical Tension 


ag 7 W000) lbs. 
an 


Thus, the horizontal component of the towline tension at 


mount A, T can be calculated from the resistance equations 


ah’ 
given by Jackson [13]. For a specific design the desired 

depth and speed of tow would be dictated by operational re- 
gGuirements. In this notional case a maximum tow velocity of 


15 knots with a design tow depth of 200 ft. will be assumed 
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as a design point. Under actual towing conditions the verti- 
cal component of the towline tension must be compensated for 
by a countering force on the control surfaces. The range of 
Bossible control surface forces would depend upon the flow 
velocity, the size of the control surface, and the angle of 
deflection. Since these conditions would be specific to each 
case, a maximum vertical tension component, eer, Of 5, 000 
lbs will be assumed with an optimum value of 7,000 lbs during 
normal operations at the design point velocity of 15 knots. 
For a given tow velocity, U, the actual velocity of the 
towed vessel and towline relative to the water, V, would be 
the algebraic sum of the towship velocity and the current 
velocity, v. It will be assumed that the local current ve- 
locity is small compared to the velocity of the towline 
through the water and that the relative velocity is equal to 
the towship velocity. It is recognized that this is a weak 
assumption for actual operations which could be conducted 
under conditions where the current velocity may be of the 
same order of magnitude as the towing velocity. However, if 
the correct relative velocity were known, it could be used 


in place of the towing velocity to improve the design. 


meee coverning Equations 

By removing an arbitrary element from the towline, as 
in Fig. 9, and establishing a local coordinate system which 
is everywhere normal and tangential to the towline axis, the 


governing equations will be: 


at/ds 


(Ws) See ouN (Oe) + F (14) 


t 
(W-B) x COS(¢) - BPs (15) 


T x d(6)/ds 


where T is the effective tension, as shown in Eq. (15a), 
which includes the hydrostatic pressure acting on the local 


cross-sectional area: 
T = TC Fae XZ eX A (15a) 


The hydrostatic component is included here to offset the 


equivalent opposing pressure term in the elemental buoyancy. 
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The effective tension, T, is the appropriate tension to be 
used in elastic elongation calculations because it accounts 
for the Poisson effect of the hydrostatic pressure. However, 
it remains an approximation because it is applied to the 
total rope volume as if it were a linear elastic solid cyl- 
inder. The working tension, Toe should be used in calculat- 
ing the specific load in the cases where strength limitations 
are of greater concern than total elongation, or for cases 
where the depth is great enough to make the pressure contri- 
bution significant. For the rope sizes and the depths of 
interest in this preliminary design the hydrostatic component 
is considered to be negligible, i.e. on the order of 0.43 
of Boy and the effective tension will be accepted as the 
working tension. This results in a more correct elongation 
model and a conservative strength model. 

Ee and F, are, respectively, the normal and tangential 


ics 
hydrodynamic drag forces, defined as: 


Bo 1. x px C xdx vee (16) 
Fy = 4 X 0. x C, Mob. <a x Vie (17) 


The equations are further simplified by defining the 


weight in water as: 


W 
W 


(W-B) = Py x d* xg x opp sere) (ns) 


The magnitude of the drag coefficients, C. and Ci 1s 
also an area of current synthetic line research and, again, 
approximations as found in Newman [14] and Springston [18] 
for the condition that flow-induced vibrations are not 


Significant, are: 


om = 1.0 (19) 


C, —tOn Os) tO, 0.05) Xx C. (20) 


3.5 = 





The mean values will be used in this preliminary design, 


Co = 10 (19a) 
n 


C, = 0.04 (20a) 


The use of thses values implies the assumption that flow- 
induced vibrations are not significant. The validity of 
this assumption will be examined after a preliminary design 
has been completed and the position of the towline rela- 
tive to the flow has been approximated. 
Where vortex shedding is present, and excites tow- 
line motions normal to free-stream velocity, the normal 
drag coefficient has been found to be as high as 3.0 [16]. 
This condition is usually found when the flow field is 
nearly uniform and normal to the axis of the test cable 
and, as such, would probably be minimal during towing opera- 
tions. The frictional drag has not been measured during such 
flow conditions, and the effect of flow-induced vibrations 
on the frictional drag coefficient is not known. 
A reference condition of near zero velocity with a 
Minimal tension of Ty) will be assumed for an elemental 
length of dsp) and diameter dy. This diameter is equivalent 
to the new rope diameter and with its new-dry breaking 
strength, Bor would be the criteria by which the rope would 
be ordered prior to incorporation into the towing system. 
When immersed in water, the nylon will "shrink", causing 
the length to decrease. Recent measurements taken from 
active U.S. Coast Guard towlines showed an average shrinkage 
Of approximately 5%. This will be the value used in the pres- 
entation where the new-dry length will be reduced by 5% to 


produce the wet reference length. Once installed into the 
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towing system and subjected to the cyclic loads, the perma- 
nent elongation, and thus the working length, of the line will 
be established. As previously stated, the permanent elonga- 
tion for wet double braid nylon is approximately 1332 at a load 
@c 20% of Ba: Further, if it is assumed that the double braid 
Seructure does not contain a significant amount of trapped air 
when loaded to T,, the assumption of incompressibility can be 
used in conjunction with the permanent elongation to provide 
the following relation between the reference length, ds,, and 


the working length, ds: 


ds, = 1.13 x ds, (21) 

A, x ds. ee ads- (22) 
& On<.5 

qd = d (a es) (23) 


The elemental working length, ds, will be used throughout 
this thesis as the basis for local elastic elongation. The 


resulting relation for elastic elongation is given by: 


e = (ds - ds_)/ds_, (24) 


ds = ds, x.1.13 x (1 + e) (25) 


0 


A second application of incompressibility combined 
with an assumption of small elastic elongation, less than 
15%, aS cited in [1] and [10], leads to the final equations 
needed to express the local diameter, d, as a function of 


S@enOriginal diameter, d,: 


Leo Ss = Aes as (26) 
W W 

Pia x qv auCle oe iy Aue d* x ds (297) 

Ze 2 
ad‘ = q. Xx ds /ds (28) 
d= 4 /(1 + e) 1/2 (29) 

where for small e: 

(1 + e)?/2 = (1 + 2/2) (30) 
eer at gc/2)) (31) 


reife te 


rs 





mee de (1/13) °°>/ (2 + e/2) | (B32 


Substitution of these relations and the relative veloc- 
ity, V, into the governing equations leads to a system of 


equations which must be satisfied by any valid model of the 


mowlLine: 
ards = We x SIN(®) + FL (33) 
T x d(%)/ds = W x COS($) - FL (34) 
Fo= 4 xp x C, x dy x Meee eee Clee / 2) x 
fey SiN (¢) x ds, x 1.13 x (1 + e/2) (25) 
es 0.5 
a ees px C. SaPi ox d x [lye 3 } /(1 + e/2)] 
mre x COS*(%) x ds, x 1.13 x (1 + e) (36) 
Mmmeeet 4 x d,* x [(1/1.13)°°*/(1 + e/2)]* x 
Soe X (p> - p.) x dsemex 1. ES ox (lame) (37) 


4.3 Method of Solution 


The solution of these equations will be achieved with 
an iterative numerical scheme using the drag at the towed 
body, the initial length of the towline, and the specified 
Minimum specific tension as the boundary conditions for a 
viable solution. A range of feasible towline angles is 
assumed at Point A, Fig. 8. The resistance of the towed 
body is either supplied as an input, or calculated, and the 
initial tension, at point A, is then provided by the rela- 


tion: 
T = R,/COS (9) (38) 


This tension is then multiplied by the minimum specific 
loading, obtained from the factor-of-safety analysis, to 


provide the desired breaking strength of the towline: 


B = T x ore (39) 
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Using regression analysis from the double braid nylon rope 
specification [10], the required breaking strength can be 


related to a nominal diameter in the form: 


doce. 241485) 2 (40) 


where Be Peewee ose awe S tne tnechnes. This 1S an approxi- 
mate relation which has a maximum error of 7%. However, the 
diameter and the breaking strength that result from the pre- 
liminary design will not coincide with standard available 
sizes; in which case, the next higher standard size should be 
used, provided that both the diameter and breaking strength 
exceed the values obtained from the design model. 

Rearranging Eq. (4) produces a relationship for the 


local elongation as a function of the local specific tension: 
eg uOe. 2 ee ee (41) 


For each assumed angle Eqs. (39) and (40) are solved for each 
element, ds, along the length of the line. The local diam- 
eter and elongation are then used as inputs for Eqs. (33) 
through (37). The change in depth as a function of elemental 
length and local angle is then calculated over the length of 


the towline by the geometric relation: 


aug) = 2(j > 1) + ds, x SIN(9) > cee gs I oe 


(j-1) 
x (1 + e) (42) 


The system of equations is iterated until 2(j) is equal to 
the desired depth of tow. A new breaking strength and diam- 
eter is calculated each time, with the final values being 
those that satisfy the initial design constraints. A Fortran 
program which can be used for preliminary sizing is shown in 
Appendix D. The results of this program for the design point 
given in Table 1 are presented in Fig. 10. 

The actual towline profile is shown in Fig. 10a. The 
profile is nearly flat due to the minimal difference between 
the densities of sea water and nylon, and due to the normal 
drag force which tends to lift the towline. The variation 


in tension along the length is shown in Fig. 10b, where it 


—39- 


can be seen that the total increase of 20,000 lbs. due to the 
hydrodynamic drag is of the same order of magnitude as the 
drag of the towed vessel. Figure 10c shows the local specific 
tension for the quasi-static model, the minimum specific ten- 
sion of 10% was an input parameter, and the resulting maximum 
specific tension is approximately 16.5%. This is slightly 
greater than the desired 15% limit given in Eq. (13a), but 
still leaves approximately 12% of the breaking strength to 
account for dynamic loading and to remain within the limits 
Smead. (13). 

Figure 10d shows the elastic elongation which has a 
local maximum of approximately 7%. When added to the 133 
permanent elongation, the total elongation of 20% is 4% less 
than the elongation at break cited in [10] and approximately 

% to 10% less than that cited by Flory [4] for cyclic load 
tests. 

The reduction in diameter due to local elongation is 
shown in Fig. 10e, and the local angle with respect to the 
Merizontal is shown in Fig. 10f. 

A typical output from the design program is given in 
Appendix E. For the selected design point and the notional 
submarine, the recommendation for a 1200-ft towline is given 
in Table 3. Using the recommended size we conclude that the 
best standard size is a 10-inch circumference, double braid 
nylon rope with dimensions also shown in Table 2. Further 


analysis of the towing system will be based on this standard 


size. 
TABLE 2 
Recommended Towline Size 
Recommended Size Standard Size 
Diameter 3.22" Diameter 3.25" 
B . 320,000 lbs. B 6227000. lbs. 
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Shap ter 
OFF-DESIGN ANALYSIS 


5.1 Low-Load Operations 


The initial selection of a design point has been the 
basis for calculating the preliminary size of the towline. 
These design point conditions are most often chosen to rep- 
resent the worst case loading conditions and thus should 
result in the selection of a rope size capable of sustaining 
both the static and dynamic loads associated with extended 
towing operations. This engineering philosophy is valid for 
most of the commonly used materials in marine engineering 
and, indeed, provides a good starting point for polymeric 
materials. However, some polymeric materials (in this case, 
nylon), exhibit vastly different behavior when the loading 
regime is varied from that chosen as a design point. A com- 
plete preliminary design must also include an examination of 
off-design performance. 

Some investigators [11] have recently found that the 
combination of relative movement between the structural ele- 
“ments of the rope and exposure to water can radically alter 
the abrasion characteristics of nylon filaments; the number 
of cycles to failure in laboratory experiments is an order 
of magnitude lower for wet-nylon yarns as compared to dry- 
nylon yarns. With comparable normal pressures and relative 
motion the same trend could be expected within the structure 
of the rope when it is immersed in water. This behavior has, 
in fact, been shown by Flory [4] in field experiments with 
eight-strand plaited nylon ropes of the size that could be 
used for towing vessels with a displacement less than 200 
tons. These abrasion effects have also been reported by 
Bitting [17] in double-braid nylon ropes used as deep water 
buoy moorings. He notes that the ropes were subjected to 
low amplitude cyclic loads for extended periods, precisely 
the condition that the lower specific load limit is meant 


to avoid. 
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It is obviously impossible to preclude exposure to 
water in a marine towing system. Thus the only parameters 
that can be varied to limit the abrasion effects are the 
material, the local normal pressure, and the extent of 
relative movement. At some point in the future nylon may be 
replaced by polyester, or some other material, but this is 
at best several years away and may not happen at all, if 
other polyester characteristics, such as elasticity and 
temperature sensitivity, negate the possible gains in abra- 
sion resistance. The remaining parameters (normal pressure 
and relative movement) can be somewhat controlled in the 
lower load regime by controlling the tension in the towline. 
However, these two parameters are not independent and cannot 
be considered separately. As the tension is reduced, the 
normal pressure between structural elements is also reduced 
and at some point the frictional forces will be too small to 
prevent the relaxation of the braided or twisted structure, 
thus allowing a greater degree of relative movement between 
the elements. During the next increasing load cycle the 
motion will be reversed and it is in this loading regime 
that abrasion resistance may have a pronounced effect. The 
lower specific load limit of 3% was imposed in the factor of 
safety analysis to account for this effect, but this limit 
only entered the design point analysis in an indirect 
manner by choosing a lower specific load limit of 10% in 
an attempt to stay above the 3% minimum during off-design 
Operations. It is during slow speed submerged or surface 
towing operations that the lower limit becomes important. 

Beene Submarine towing scenario used here the effects 
of varying the tow velocity are depicted in Figs. lla 
through lle. The initial tension in the towline at the 
connection point on the submarine is a function of velocity, 
which controls the hydrodynamic drag of the submarine, and 
depth of tow which dictates the initial angle of the tow- 
line. If, for the moment, we ignore the physical limitations 
of depth of tow and the vertical component of the towline 


tension, the depth of the submarine, in theory, can be 
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adjusted to create an angle which maintains a minimum tT of 

3% in the towline which may inhibit the internal abrasion. 
The resistance of the notional submarine as a function of 
velocity is shown in Fig. lla, and the tow depth that would 
maintain the required minimum specific tension is shown in 
Fig. llb. The resulting initial towline angle is shown in 
Fig. lic. As expected, the angle approaches a maximum of 90° 
as the velocity approaches zero; at this point, the depth of 
tow is equal to the stretched length of the towline. Figure 
lid depicts the variation of t along the towline, and it can 
be seen that the maximum value at the surface tow point is 
well below the limit given in the factor of safety analysis. 
However, at this point we must invoke reality and examine 

the vertical component of the towline tension which must be 
countered by the control surfaces of the submarine. Since 
the forces acting on the control surfaces are also a function 
of flow velocity, there will be a physical minimum speed and 
maximum depth of operation for each case. If, for the pur- 
pose of illustration, we assume that the control surface 


force is a V-squared function of the form: 


= 2 
Fae = Ce V (43) 


For the notional submarine C is 67 and the control surface 
limit is reached between 7 and 8 knots, as shown in Fig. lle. 
This translates into a depth limit of approximately 400 feet 
mir ig. lib. 

The depth of tow is further limited by the maximum 
Operational depth of the specific submarine being towed; the 
designer would thus have to impose this limit of Fig. 11b. 
The resulting maximum depth; either the operational depth 
limit, or the control surface limit, would then reflect the 
Minimum velocity that would maintain the required minimum 
Specific loading. 

The decrease in the abrasion resistance of nylon when 
exposed to water has been well documented. However, it should 


be emphasized that the exact magnitude of the lower limit for 
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a particular rope structure is not yet established. Contin- 
uing research into the wet abrasion characteristics of nylon 
may either raise or lower the limit in relation to the 3% of 
breaking strength cited by Bitting and used here. In general 
terms, it can be stated that any increase in the limit would 
only further restrict the operating envelope and any slight 
decrease would offer only minimal gains because the limit is 


already very low. 


meee ocurface Towing 
It is highly unlikely that all submarine towing opera- 


tions would be conducted with the submarine below the surface 
of the water. While entering and leaving port, or during 
emergency situations, the submarine would most probably be 

on the surface. Analysis of this scenario requires a knowl- 
edge of the positioning of the towing hardware on the bow of 
the submarine. If the tow point 1s above the water, the scope 
is short, and the towline is kept taut, the hydrodynamic drag 
of the towline would not be a factor, and the initial towline 
angle would be very nearly zero. In this case, the minimum 
tension is equal to the drag of the submarine and would be 
below the 3% minimum when the drag of the towed vessel is 

less than 3% of the breaking strength of the towline. For 

the notional submarine and the 3.25"-diameter towline, this 
would occur at a tow speed of approximately 7 knots. tThis 

is Significant since all towing operations in restricted 
waters and during many emergency conditions would be conducted 
at less than 7 knots. 

For the case where the tow point is below the surface of 
the water, a depth of five feet is assumed here, hydrodynamic 
drag on the towline must be accounted for. Since the towline 
is nearly parallel to the velocity vector, only the tangential 
drag component will be significant. Figures Al through A4 
in Appendix A, present the plots for towline profiles, tow- 
line angle, specific tensions, and elastic elongations at 
various towing velocities for surface towing with a submerged 


tow point. At 3 knots, Fig. A3 shows that the specific 
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tension 1s everywhere less than 3% and Fig. A4 shows that the 
elastic elongation is less than 2%. This is precisely the 
condition that would exist in most cases for entering and 
leaving port and it is a condition where abrasion could be sig- 
nificant. The same characteristics are shown for a tow ve- 
mocaty of 6 knots. The initial t of 2.8% is just below the 
recommended minimum and depicts an acceptable towing condition. 
At 9 knots t has increased to an average value of 10% with an 
associated average elastic elongation of 5.7%; both of these 
values are acceptable, but they are now approaching the maxi- 
mum limit for total elongation, 24%, where the elastic elon- 
gation at the towing vessel is added to the assumed permanent 
elongation of 13%. The 12-knot curve is above the upper limit 
for t (Eq. (1l3a))}, and represents the maximum surface towing 


velocity for calm water where dynamic loads are not signifi- 


eant. 
TABLE 3 
Notional Surface Resistances 
Velocity in knots Resistance in pounds 
3 4,000 
6 9,000 
9 30,000 
12 60,000 
15 7s, 008 


The surface resistance characteristics given in Table 3 
have been assumed for a notional submarine; an actual towing 
analysis would be based on known resistance data, or scale 
model tests. The notional submarine used here simply illus- 
trates the possibility that a towline sized for a submerged 
velocity of 15 knots may be more severely limited for surface 
towing conditions. Here, we have an approximate surface 
Maximum of 10 knots which could decrease even further if the 
assumed dynamic amplification factor was not large enough to 
account for wave excitation of both the towing and towed 


vessels. 
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5.3 Flow-Induced Oscillations 


Flow-induced oscillations of long cylinders in a fluid 
medium is a dynamic phenomenon which has been investigated 
by, among others, Sarpkaya [15], for a stationary bluff body, 
and by Vandiver [16] for lcng flexible cylinders. Considera- 
tion in a quaSi-static analysis 1s only to the extent that 
such vibrations have been reported to increase the normal 
drag coefficient by as much as 200% from 1.0 to a maximum of 
3.0. Such increases have been measured in laboratory and 
field experiments where the flow was normal to the axis of the 


cylinder and under the following conditions: 


a a 1 
pane Spatially uniform cylinders 
Pie opat load ly Unttenm flows 


A Low modal density. 


y 


Under these conditions it is possible to approximate the 
normal drag coefficient with some degree of confidence by 
using an empirical relation proposed by Griffin and modified 


by Vandiver in the form: 


= 0,65 
CC, = Cue * [1+1.043 (2¥_ ./d) } (44) 
where oa Heethe £igid ¢cylinder normal drag coefficient 

is the RMS ampli- 
rms 


tude of cable vibration. If, for a worst case, we assume 


that has been used to this point and Y 


that the towline can be approximated by the above conditions 
mmoamthat the amplitude of vibration is 0.7 x A, similar to 
that reported for the pipe by Vandiver [16], we find that the 
Co maximum 1S 2.8. A normal drag coefficient of this mag- 
nitude would obviously result in a substantial increase in 
the size of the towline which would be capable of sustaining 
the maximum loads at the end connected to the towing vessel. 
Since the resistance of the towed vessel would remain the 
same for a given speed, it would also make it much harder to 
remain above the lower load limit at that end of the towline. 
The Reynolds number, Rye for the actual towline is shown 


in Fig. 12a. It has been calculated, based on the local 
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mormal velocity, Fig. 12b and the local diameter, Fig. 12c. 
From these figures we see that the first and second con- 
ditions, R, <= 10° and a spatially uniform cylinder, have 
met, Dut that the third condition, spatially uniform flow, 
is not met. Since the flow field is not uniform, we expect 
a very short vortex correlation length with a random vibra- 
tion response [16]. Under these conditions, the amplitude 
of vibration and the associated drag are substantially re- 
duced. CA approaches Cao and a magnitude of 1.0 to 1.2 is 
found to be a good approximation for sizing the towline. 

Figures Bl through B6, showing the effects on the 
required size and the equilibrium condition due to increas- 
ing the normal drag coefficient from that assumed for design, 
i.e. to a maximum of 2.5, are presented in Appendix B for the 
design point of 15 knots and 200 feet. In practice one would 
probably purchase the standard commercial size above that 
recommended by the design. A design based on a Co os bie! ey 
would result in a 10-inch circumference towline, as was al- 
ready recommended using a Cc. SteiyO.e im Contrast, a design 
based ona C. of 2.5 would result in the selection of a 1l2- 
inch circumference towline. This variation is not minor, 
Since the increase in Be is of the order of 40% and it would 
become impossible to keep the specific load above the lower 
limit for slow speed tows or during any configuration where 
flow-induced vibrations decreased or stopped and the normal 
drag coefficient decreased to approximately 1.0. 

For all cases, the tangential drag coefficient was 


assumed to be constant at 0.04. 


fee Dry- Versus Wet-Nylon Function 


Barlier we examined the elastic elongation characteris- 
tics of double braid nylon and reached the conclusion that 
under high cycle load conditions the DNF better represented 
the data presented in the literature. It was also stated 
that the WNF would be more appropriate for a new line being 


placed in a towing system, and that it would provide an upper 
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limit on in-service elongation. However, the point at which 
the design switches from the WNF to the DNF has not been 
established. From Figs. 2b and 3 it would appear that this 
change could be made as early as the first few hundred 
cycles, but cyclic elongation data for double braid nylon 
rope used in marine towing systems is very limited, and a 
conservative approach would be to use the WNF for the first 
few thousand cycles. This translates into 4 to 8 hours of 
towing under normal operating conditions. The implications 
of using the WNF are presented in Appendix C, for operations 
at or near the 15-knot and 200-foot design point, and for 
Slow speed, deeper operations at 3 knots and 600 feet. 

As might be expected, the towline profile, towline ten- 
Sion, and specific tension are practically the same for both 
the DNF and the WNF, Figs. Cl through C3, and Cll through 
C13. In contrast, the elastic elongation in Figs. C4 and Cl4 
show a Significant increase when the WNF is used. Because of 
the increase in elongation there is a concurrent decrease in 
local diameter as shown in Figs. C5 and C15. At the design 
point, the maximum elastic elongation is approximately 133, 
which, when combined with the assumed permanent elongation 
of 13%, implies a total elongation of 26%. This is an un- 
acceptable value which leaves little or no elongation margin 
for dynamic loads. It indicates that a new towline has an 
operating point substantially below that of a used, but un- 
damaged, towline. 

The remaining figures show only minor effects from the 
increased elongation, and it can be concluded that upper 
range of the operating envelope is limited by local elonga- 
tion rather than strength, as reflected in the specific 
loading and thus the factor of safety. However, it can also 
be concluded from these figures that a towline size based on 
the WNF, which would be substantially larger, would be very 
prone to the effects of internal abrasion, since it would 
continually operate at much smaller specific loads. The 
specific loading conditions shown in Fig. Cl3 would be typical 


of such operations even at higher towing speeds. 
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Chapece Ve 
CONCLUSIONS 


mee Controlling Material Characteristic 


In the preceding discussion it has become apparent that 
recent research [11] which has revealed the susceptibility 
of wet nylon to abrasion and the apparent greater dominance 
of this characteristic under low load conditions Piy |. asethe 
most severe constraint on designing a nylon towline to oper- 
ate over a wide range of loading conditions in the marine 
environment. It is theoretically possible to increase the 
meee Or a towline to account for larger static or dynamic 
loads, or to account for a greater degree of uncertainty, 
but, because of internal abrasion that has been observed in 
nylon lines, such a design approach simply changes the 
mechanism by which the towline fails and does not decrease 
the likelihood of failure. 


Meee controlling Structural Characteristic 


The controlling structural characteristic in this de- 
Sign model is the assumed form of the elastic elongation 
function. The DNF and the WNF result in significantly dif- 
ferent elongation behavior. Because of established abrasion 
characteristics of wet nylon, neither function can be used to 
produce a conservative design over the total projected loading 
regime of a marine towline. Both of these functions are 
based on a minimal number of cycles and even though they ap- 
pear to bound cyclic load data presented by other authors 
(for what are assumed to be wet testing conditions), it is 
noted that there are no truly wet tests where the rope has 
been continuously immersed in water, available for compari- 
son. As a result, the applicability of the elastic elonga- 
tion function to very high cycle wet loading conditions is 
questionable. Future research may provide better models 
which result in a much less constrained design. Until this 


research has been completed, some form of the elongation 
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behavior presented here must be used in conjunction with an 
experienced seaman's knowledge of what has worked well in the 


past. 


SeoeecoOntrolling Loading Parameter 


A review of Figs. C2 and C10 shows that the appreciable 
increase in tension along the towline is primarily due to the 
frictional (tangential) drag of the towline. The magnitude 
of F, is a reflection of the tangential drag coefficient, C,- 


t 
The value of C, used here was suggested by Triantafyllou and 


supported by as alee [15], and is probably very close to 
the actual value, but it appears that there are very little 
data specifically for the frictional characteristics of syn- 
thetic ropes. It would require only relatively small 
variations of C, EOmOEOGdNce ta falmlyusiqniticant change in 
the maximum local loading and the associated maximum local 
elongation. This change would be beneficial if C, were found 
to be less than that used here and would be detrimental if it 
were found to be greater. In either case, it would substan- 
tially improve the confidence in the final selection of the 


rope size. 


6.4 Summary 


One possible method of selecting the size of a double 
braid nylon rope for a marine towline has been presented 
here. This method is primarily based on behavior character- 
istics presented in the literature. It was not possible to 
reach a definitive conclusion because the model was event- 
ually controlled by assumed values for the minimum acceptable 
ieaaing COndition and the appropriate hydrodynamic frictional 
drag coefficient. These assumed values have a basis in the 
literature, but further research is needed to establish their 
exact magnitudes for synthetic ropes of various materials 


Pmoy Or Constructions. 
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APPENDIX D 


ml. Design snaorengreuin 


CI" CQ Cy Ca" Cl ar Gar Ca Ca cle G2) 


aqgaaaana 


aaa 


gqagaa 


THIS PROGRAM JS USED TO PREDICT THE SIZB OF A DOUBLE BRAID 
NYLON ROPE WHICH CAN BE USED AS A MARINE TOWLINE 

THE ALGORITHM IS BASED ON A SUBMARINE TOWING SYSTEM WHERE 
THE DESIGN POINT IS FOR A SUBMERGED TOWING CONDITION AT A 
PREDETBRMINED DEPTH AND SPEED 


INPUT PARAMETERS CONSIST OF THE GEOMETRY OF THE TOWED VESSEL, THE 
LENGTH OF THE DESIRED TOWLINE AND THE DESIRED MINIMUM LOAD AS A 
PERCENTAGE OF THE RATED BREAKING STRENGTH, A MINIMUM VALUE OF 0.08 
TO 01 IS STRONGLY RECOMMENDED TO AVOID POOR OFF-DESIGN BEHAVIOR 


IMPLICIT RBAL(L-N) 
DIMENSION (1200), T(1200), PHI(1200),2(1200),X(1200), Y(250),ST(25 
0), DI(1200), TH(1200) 


THE VALUB CHOSEN FOR CD HERE IS 10, THIS PARAMETER HAS A 
SIGNIFICANT EFFECT ON THE UPPER LOAD LIMIT AND SHOULD BE CHOSEN 
WITH CARE 


DATA RHOW RHOLPI,CD/20,2 209,3.141593,1 0/ 
WRITE(?,101) 
RBAD(*,102) L 
WRITE(?,103) 
RBAD(*,102) D 
WRITE(*,105) 
RBAD(*,102) VK 
WRITE(*,104) 
RBAD(* 102) TD 
WRITE(*,106) 
RBAD(*,107) IL 
WRITE(*,108) 
RBAD(*,102) SLMIN 


THE FOLLOWING LINE ASSUME A CONSTANT SHRINKAGE OF 54% FOR NYLON 


SHRK=FLOAT(IL)- 05*FLOAT(IL) 
CT=004*CD 
IL=INT(SHRK) 


Ai THRU RT PROVIDE AN APPROXIMATION OF THE INITIAL TENSION IN 
THE TOWLINEB FROM THE RESISTANCE OF THE SUBMARINE 


Al=(VK*L*1 69)/.0000129 
A?=(ALOGIO(A!)-2)**2 
A3=L/D-1 3606 

A4=075/A2+ 0008+ 000789/A3 
AS=PI*D**#2*A3*A4 


SiGe 
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A6=A5+9 + 001*L'D 
DHP=.00872*VK**3*A6 
RT=(DHP*550)/(1 6886*VK) 
WRITE(*,200) RT 


THE INITIAL ANGLES ARE SET BQUAL TO THE MAXIMUM PHYSICALLY 
POSSIBLE FOR THIS TOWING GEOMETRY 


TAU1=SLMIN 
UPHI=PI/2 -0! 
LPHI=.0! 

DS=}. 
MPHI=(UPHI+LPHI)/2 
TAU=TAU! 


INITIAL BREAKING STRENGTH AND DIA BASED ON SAMSON DATA 


BSL=RT/(SLMIN*COS(MPHI)) 
DIL=(BSL/34148 5)**.5258 


HYDROSTATIC COMPONENT OF TENSION 


TH(1)=.3491*TD*DIL**2 
DO 10 I=1,3 


INITIAL EFFECTIVE TENSION AS A FCN OF ANGLE 


IF(I~2) 51,5253 
T(1)=RT/COS(LPHI) 
PHI(1)=LPHI 

GO TO 55 
T(1)=RT/COS(UPHI) 
PHI(1)=UPHI 

GO TO 55 
T(1}=RT/COS(MPHI) 
PHI(1)=MPHI 


INITIAL ELASTIC ELONGATION FROM DNF 
B(1)=.2119*(TAU1)*# 5848 

DI(1}=DIL*.941 /(1+B(1)/2) 

Z(1}=0 

X(1)=0 

TE=E(1) 

BEGIN ITERITIVE SOLN OF GOVN EQN 


DO 20 J=2IL 


ai = 





aa 


20 


99 
61 


62 


ESTABLISH SIMPSONS MULTIPLIERS FOR TOTAL BLONGATION 


IF((2*INT(J/2)) LT.J) TRAC=2 
IF((2*INT(J/2)) BQ J) TFAC=4 
IF(J EQIL) TFAC=! 

IN=! 

LOCAL ELAS BLONG 
B(J}=.2119°(TAU)** 5848 
TE=TE+E(J)*TFAC 

LOCAL DIA. 

DI(J)=DIL* 941/(1+B(J)/2) 


LOCAL WEIGHT IN WATER 
WB=PI/4 *(DI(J)/12 )**2*322*(1+B(J))*(RHOL-RHOW)*! 13 
LOCAL NORMAL DRAG FORCE 
FD1=5*RHOW*CD*(DI(J)/12 }*(WK*1 6886)**2 
FD2=SIN(PHI(J-IN))**2*(1+B(J))*113 

FD=FD1*FD?2 

LOCAL TANGENTIAL DRAG FORCE 
FT1=5*RHOW*CT*PI*(Di(J)/12 )*(VK*1 6886)**2 
FT2=COS(PHI(J-IN))**2*(1+B(J))*113 

PT=FT1*FT2 

LOCAL ANGLE 
PHI(J)=PHI(J-IN)+DS/T(J-IN)*(-FD+WB*COS(PHI(J-IN))) 
DEPTH OF ELEMENT 
2(J)=Z(J-IN)+DS*SIN(PHI(J-IN))*(1+B(J))*1.13 

LOCAL HYDROSTATIC TENSION 

TH(J}=0 3491*(TD-2(J))*DI(J)**2 

LOCAL EFFECTIVE TENSION 
T(J)}=T(J-IN)+DS*(FT+WB*SIN(PHI(J-IN))) 
TAU=T(J)/BSL 
X(J)=X(J-IN)+DS*COS(PHI(J-IN))*(1+B(J))*113 

IM=J 

CHECK FOR IMPOSSIBLE SOLUTIONS 
IF((Z(J)-TD).GT.10) GO TO 59 

IF((TD-Z(J)). GT (2*TD)) GO TO 59 

CONTINUE 

CHECK FOR SOLUTION BASED ON FINAL DEPTH AT STRETCHED LENGTH 
IF(1-2) 61,6263 

SZL=2(JM) 

SXL=X(JM) 

FZL=TD-SZL 

BRROR CONDITION ON LOWER ANGLE 

IF(FZLLT 0) GO TO 1000 

GO TO 10 

SZU=Z(JM) 

SXU=X(JM) 

FZU=TD-SZU 

ERROR CONDITION ON UPPER ANGLE 

IF(FZU GT 0) GO TO 1010 
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30 


80 


260 


101 
102 
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104 
105 
106 


GO TO 10 
SZM=Z(JM) 

SXM=X(JM) 

FZM=TD-SZM 

CHECK FOR DEPTH SOLN CLOSE TO SURFACE 

IF(ABS(FZM)LT.1) GO TO 90 

CONTINUE 

CHECK TO SEE IF LIMITS ARE ON OPPOSITE SIDES OF SOLN 

IF NOT WRITE ERROR MESSAGE 

IF((FZU*FZL)GT0.) GO TO 80 

DPHI=UPHI-LPHI 

ERROR CHECK ON ANGLES 

IF (DPHIEQ.0) GO TO 95 

CHECK FOR SOLN BASED ON CONVERGING ANGLES 

IF(DPHILT 0001) GO TO 90 

MODIFY UPPER OR LOWER ANGLE LIMIT FOR NEXT ITERATION 
IF((FZU*FZM)GT.0) UPHI=MPHI 


IF((FZL*FZM)GT.0) LPHI=MPH! 

SCREEN WRITE TO MONITOR CONVERGENT BEHAVIOR 
WRITE(*,260) FZL,FZM FZU LPH], MPHIUPHI 

START NEXT ITERATION IF NEEDED 

GO TO 15 

SLMAX=T(JM)/BSL 

TE=(DS/3 *TB+13*FLOAT(IL))/FLOAT(IL) 

WRITE FINAL SOLN TO DATA FILE 
OPEN(16,FILE="TOWNYLON DAT’ STATUS="NEW’) 
WRITE(16,225) VK. RT.TDIL,TR 

WRITE(16,220) DIL,BSL,SLMINSSLMAX 

WRITE(16,230) 

1V=5 

DO 30 K=1 ILIV 

IF(T(K)LE9.) GOTO 30 

Y(K)=2(K)-TD 

ST(K)=T(K)/BSL 

WRITE(16,240) B(K),T(K), TH(K),ST(K),DI(K),PHI(K), Y(K),X(K) 
CONTINUE 

CLOSE(16) 

GO TO 95 

WRITE(?,210) 

WRITE(*,260) FZL,FZM FZU,LPHI,MPHI,UPHI 

FORMAT( °? FZL B92? FZM °,B92? FZU B92’ LPH] °,B114, 
MPHI ?.B114” UPHI ’,B114/) 

FORMAT( > INPUT LENGTH OF TOWED VESSEL °) 
FORMAT( F62) 

FORMAT(’ INPUT DIAMETER OF TOWED VESSEL °) 
FORMAT(’ INPUT DESIGN DEPTH OF TOW °) 
FORMAT(’ INPUT DESIRED DESIGN VELOCITY IN KNTS °”) 
FORMAT( ° INPUT LENGTH OF TOWLINE TO THE NEAREST 10 FT,AS AN INT 


1EGER °) 


270= 
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108 


200 


210 


220 
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230 


240 


290 


1000 
270 


1010 
280 
95 


FORMAT( 16) 

FORMAT( ’ INPUT DESIRED MINIMUM SPECIFIC LOADING AT TOWED VESSEL 
©) 

FORMAT( ’ RESISTANCE IN LBF 2, F102) 

FORMAT( ’ ANGLE NOT WITHIN RANGE OF INITIAL ESTIMATE") 

FORMAT( ’ ROPE DIAMETER IN INCHES > F64’ NEW-DRY BREAKING STR 


1ENGTH LBS ?,F102/ MINIMUM SPECIFIC LOAD °F63” MAXIMUM SPEC 
2IFIC LOAD °,F63) 


FORMAT(’? TOW VELOCITY °*F61,’ RESISTANCE ?F81? TOW DEPTH 


1° F61/) TOWLINE LENGTH 716? TOTAL ELONGATION °.F83) 


FORMAT(’ ELONGATION TENSION THYDRO TAU DIA A 


INGLE DEPTH REACH’) 


FORMAT( ?F65,4X,F103,1X,F103,2X F65,4X,F643X,F64 2X F8 


1.3,3X,F9 3) 


FORMAT( ’ ?,2F62) 

GO TO 95 

WRITE(* 270) FZL 

FORMAT( ’ LOWER LIMIT IS ABOVE WATERLINE! ’.B146) 
GO TO 95 

WRITE(*,280) FZU 

FORMAT( ’ UPPPER LIMIT IS BELOW WATERLINE! ’.B146) 
END 


D2. Analysis Program 


OO A? BP OS BOO eo NOS Wie ee ee 


THIS PROGRAM IS USED TO ANALYSE THE OFF-DESIGN PERFORMANCE 
OF A PRE-SIZED NYLON TOWLINE, IT CALCULATES BOTH THE DNF AND 
WNF LIMITS FOR HIGH CYCLE AND NEW BEHAVIOR RESPECTIVELY 


THERE ARE TWO DATA FILES, ONE FOR DNF BEHAVIOR AND ONE FOR 
WNF BEHAVIOR; GENERAL DATA IS INCLUDED AS A HEADER IN THE 
WNF OUTPUT FILE 


THE DRAG OF THE VESSEL BEING TOWED MUST BE KNOWN AND PROVIDED 
AS INPUT TO THE PROGRAM 


THE OUTPUT INCLUDES THE LOCAL REYNOLDS NUMBER, NORMAL AND 
TANGENTIAL VELOCITIES FOR LOCAL FLOW CHARACTERISTICS 


IMPLICIT RBAL(L-M) 

DIMENSION E(2,1200), T(2,1200), PHI(2,1200),2(2,1200), X(2,1200) 
DIMENSION Y(2,250),ST(2,250),DI(2,1200), SLMAX(2),TB(2) 
DIMENSION FD(2,1200) FT(2,1200), VN(2,1200), VT(2,1200) RN(2,1200) 
DIMENSION TH(2,1200) 

DATA RHOW,RHOL,PI/2.0,2.209,3.141593/ 

WRITE(* 101) 


Soi 
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of 
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95 


RBAD(*,102) DIL 
WRITE(*,103) 
RBAD(*,102) BSL 
WRITE(*,105) 
READ(*,102) VK 
WRITE(*,104) 
RBAD(*,102) TD 
WRITE(*,106) 
RBAD(*,107) IL 
WRITE(*,198) 
RBAD(*,102) RT 
WRITE(*,109) 
RBAD(*,102) CD 
CT=0.04*CD 
SHRK=FLOAT(IL)-.0S*FLOAT(IL) 
IL=INT(SHRK) 
FD(1,1)=0 

FD(2,1)=0 

FT(1,1)=0 

FT(2,1)=0 

DO 25 |J=1,2 
UPHI=PI/2.-01 
LPHI=.01-P1/2 

DS=1. 
TH(1,1)=.3491*TD*DIL**2 
2 1)=T H(1,1) 
MPH]=(UPHI+LPHI)/2 
DO 10 1=1,3 

IF(I-2) 51,52,53 
T(1J,1L)}=RT/COS(LPHD 
Pail) 1)=LPHI 

GO TO 55 
T(1J5,1)=RT/COS(UPHI) 


PHI(1J,1)=UPHI 

GO TO 55 

T(1J,1)=RT/COS(MPHI) 
PHI(IJ.1)=MPHI 

TAU1=T(1J,1)/BSL 

TAU=TAU! 

IF(IJBQ1) B(1J,1}= 307*(TAUL)** 518 
IF(1J.BQ.2) BIJ1)=2119*(TAUL)** 5848 
DI(1J,1)=DIL*.941 /(1+B(1J,1)/2) 
Z(1J,1)=0. 

X(1J,1)=0 

TE(1J)=B(1J,1) 

DO 20 J=2IL 

IF((2*INT(J/2)) LT.J) TPAC=2 
IF((2*INT(J/2)) BQJ) TFAC=4 
IF(JEQIL) TFAC=1. 


=O 





20 
99 
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62 


63 


10 


90 


IN=! 

IF(1J.BQ1) B(1J,J)= 307*(TAU)** 518 

IF(1J BQ 2) B(IJ,J)= 2119*(TAU)** 5848 
TE(IJ)=TE(1J)+E(1,J)*TPAC 
DI(1J,J)=DIL* 941 /(1+B(1J,J)/2) 

WB=PI/4 *(DI(IJ,J)/12 )**2*32 2*(1+B(1J,J))*(RHOL-RHOW)?1 13 
FDi=5*RHOW*CD*(DI(1J,J)/12.)*(WK*1 6886)**2 
FD2=SIN(PHI(IJ,J-IN))**2*(1+B(1J,J))*113 
FD(1J,J)=FD1*FD2 
FT1=5*RHOW*CT*PI*(DI(1J,J)/12.)*(VK*1 6886)**2 
FT2=COS(PHI(IJ,J-IN))**2*(1+B(1J,J))*113 
PT(J)=FPT1°FT2 
PHI(1J,J)=PHI(1J,J-IN)+DS/T(IJ,J-IN)*(-FD(1J,J)+WB* 


1COS(PHI(1J,J-IN))) 


Z(1J,J)=Z(1J ,J-IN)+DS*SIN(PHI(IJ,J-IN))*(1+B(1J,J))*1.13 
TH(1J,J)=.$491*(TD-Z(1J,J))*DI(1J,J)**2 
T(1),JJ=T(1J,J-IN)+DS*(FT(1J,J)+WB*SIN(PHI(IJ,J-IN))) 
TAU=T(1J,J)/BSL 
X(1J,J)=X(1J,J-IN)+DS*COS(PHI(1J,J-IN))*(1+B(1J,J))*1.13 
IM=J 

IF((Z(1J,J)-TD)GT.10.) GO TO 59 
IF((TD-Z(1J,J))GT(2*TD)) GO TO 59 

CONTINUE 

IF(I-2) 61,62,63 

SZL=2(1J,JM) 

SXL=X(IJ,JM) 

FZL=TD-SZL 

IF(FZLLT.0) GO TO 1000 

GO TO 10 

SZU=2(1J, JM) 

SXU=X(IJ,JM) 

FZU=TD-SZU 

IF(FZU GT0) GO TO 1010 

GO TO 10 

SZM=2(1J JM) 

SXM=X(IJ,JM) 

FZM=TD-SZM 

IF(ABS(FZM)LT.01) GO TO 90 

CONTINUE 


IF((FZU*FZL)GT 0) GO TO 80 
DPHI=UPHI-LPHI 

IF (DPHIEQ0) GO TO 95 

IF(DPHLLT 0.0001) GO TO 90 

IF((FZU*FZM) GT 0.) UPHI=MPHI 
IF((FZL*FZM)GT0) LPHI=MPHI 
WRITE(*,260) FZL,FZM,FZU,LPHI,MPHI,UPHI 
GO TO 15 

SLMAX(IJ)=T(1J,JM)/BSL 
SLMIN=T(1J,1)/BSL 


aoe 





TE(1J)=(DS/3 *TE(1J)+ 13*FLOAT(IL))/FLOAT(IL) 
25 CONTINUE 
DO 35 Il=1,2 
DO 45 JJ=1IL 
VN(II,JJ)=VK*SIN(PHI(II,JJ)) 
VT(I,JJ)=WK*COS(PHI(II,JJ)) 
RN(II,JJ)=(WN(II,JJ)*1 6889*DI(I1,JJ))/(1 47B-5*12) 
45 CONTINUE 
35 CONTINUE 
OPEN(16, FILB=’"TOWANW DAT’ STATUS=’NEW’) 
OPEN(17, FILE=’TOWAND DAT’ STATUS=’NEW’) 
OPEN(18,PILE=’FVRW DAT’ STATUS=’NEW’) 
OPEN(19 FILE=’F VRD DAT’. STATUS=’NBW’) 
WRITE(16,225) VK,RT,TDJL,TE(1),TE(2),CD 
WRITE(16,220) DIL,BSLSLMIN, SLMAX(1),SLMAX(2) 
WRITE(18,225) VK RT, TDIL,TE(1),TE(2), CD 
WRITE(18,220)DIL,BSL,SLMIN SLMAX(1), SLMAX(2) 
WRITE(18,232) 
WRITE(19,238) 
WRITE(16,230) 
WRITE(17,235) 
IV=10 
DO 30 K=1L ILIV 
Y(1,K)=2(1,K)-TD 
Y(2K)=2(2,K)-TD 
ST(1,K)=T(1,K)/BSL 
ST(2,K)=T(2,K)/BSL 
WRITE(16,240) B(1,K),T(1,K), TH(1,K),ST(1,K),DI(1,K),PHI(1,K),Y( 
11K), X(LK) 
WRITE(18,242) VN(1K), VT(1,K),RN(1,K),FD(1,K),FT(K) 
WRITE(17,245) E(2,K),T(2,K), PH(2,K),ST(2,K),DI(2,K),PHI(2,K), Y( 
12,K),X(2,K) 
WRITE(19,242) WN(2,K), VT(2,K), RN(2,K),FD(2,K),FT(2,K) 
30 CONTINUE 
CLOSE(16) 
GO TO 95 
80 WRITE(* 210) 
WRITE(*,260) FZL,FZM,FZU,LPHI,MPHI,UPHI 
260 FORMAT(? FZL B92? FZM ’,B92> FZU °>B92? LPHI B11 4? 
1 MPHI >Bi14? UPHI B114/) 
101. FORMAT(’ INPUT DIAMETER OF TOWLINE °) 
102 FORMAT( F103) 
103. FORMAT(? INPUT NDBS OF TOWLINE °) 
104 FORMAT(° INPUT DESIGN DEPTH OF TOW ”) 
105  FORMAT(° INPUT DESIRED DESIGN VELOCITY IN KNTS ”) 


106 FORMAT(’ INPUT LENGTH OF TOWLINE TO THE NEAREST 10 FT,AS AN INT 
1EGER ”) 

107 FORMATT( 16) 

108 FORMAT(’ INPUT RESISTANCE OF TOWED VESSEL °) 


Oe 





109 
210 
220 


225 


230 


232 
235 


238 
240 


242 
245 


1000 
270 


1010 
280 
95 


FORMAT(’ INPUT DESIRED DRAG COBFFICIBNT, CD ’) 
FORMAT( ’ ANGLE NOT WITHIN RANGE OF INITIAL BSTIMATB?) 
FORMAT(’ ROPE DIA IN INCHES °.F64” NBW-DRY BREAKING STRENGTH 
1 LBS 2, F102/ MINIMUM SPECIFIC LOAD °?,F63/? WE MAXIMUM SPECIF 
21C LOAD ?,F63, DBE MAXIMUM SPECIFIC LOAD ’F63) 
FORMAT(’ TOW VELOCITY . ’\F61, RESISTANCE ?,F8!’ TOW DEPTH 
1° F6.1/? TOWLINE LENGTH °J6” WE TOTAL BLONGATION ?F83,’ DE 
2TOTAL ELONGATION °.F83” NORMAL DRAG COEFFICIENT °,F8 4) 
FORMAT(’ WE BELONG WETEN WETH WETAU WEDIA WEA 
INGLE WE DEPTH WE RBACH’) 
FORMAT(’??? VNW 3 VTW RN FDW FTW ’) 
FORMAT(’ DE BELONG DETEN WBETH DBTAU DEDIA DEA 
INGLE DB DEPTH DE RBACH’) 


FORMAT(’?? VND  VTW RN FDD FTD ’) 

FORMAT( ?’,F65,4X,F103,1X,F8.3,3X,F65,4X,F6 4.3X,F6 4,2X,F8 
13, 4X F9 3) 

FORMAT( '’F63,2X.F6 3,2X,B103,1XF7.31X,F73) 

FORMAT(? F65.4X,F103,1X.F83,3X F65,.4XF6 4,3X,F6 4,2XF8 
13, 4X,F9 3) 

GO TO 95 

WRITE(*,270) FZL 

FORMAT( ’ LOWER LIMIT IS ABOVE WATERLINE! ’,B146) 

GO TO 95 

WRITE(*,280) FZU 

FORMAT(’ UPPPER LIMIT IS BELOW WATERLINE! ’,B146) 


END 


OA vas 





ANALYSIS PROGRAM OUTPUT 


TYPE A8122915, DAT 


TOY VELOCITY : 
TOWLINE LENGTH : 


15.8 RESISTANCE : 
1149 WE TOTAL ELONGATION ° 


97 BORMAL BRAG COEFFICIENT ° 
ROPE DIA IN INCHES °3.2586 NEN-DRY BREAKING STRENGTH LES : 322899, #9 
MINIMUM SPECIFIC LOAD > 9.899 
WE MAXIMUM SPECIFIC LOAD : 9.183 DE MAXIMUM SPECIFIC LOAD : 9.182 


HE ELONG 


28725! 
99262 
moe? 
BE Oo2 


WW ae 


0297357 
09492 
89437 
39595 
89541 
£9579 
89689 
59643 
89678 
SEL 
89746 
99779 
£9913 
89947 
. 29231 
29914 
09948 
39981 
12614 
19947 
{9989 
19113 
18146 
18179 
» 19212 
19285 
19277 
239313 
» 19342 
18375 
19497 
2 19439 
13473 
» 18583 
183 
Biecs? 


WE TEN 
31646, 256 
31979, 342 
521t2.632 
$2346,274 
32009, 298 
$2914, 148 
33943, 441 
SERS ee 


hats hey [ata 
WO ee doe 


3378/.648 
34222, 941 
34458. 422 
54694. 982 
34929914 
$0169. 918 
30482. 982 
$3458. 496 
35874, 987 
Sole! 0 
96543.293 
Segdderd 
sedaz.c/7 
37859.477 
97296.913 
Sravaniecs 
37771.871 
52999. 586 
39247. 422 
ce Dogue 
39722, 153 


awe tus 


38961 .641 
aie ate 
57438. 544 
59676.859 
39915. 489 
49154. 283 
49393, 831 
4$631.957 
49378, 922 
411g. 
41349.5 


ro * 


4 


= 
i 
| 


aE TAU 


99828 
89999 
99973 
19845 
Letts 


eervauy 


19848 
. 19921 
18994 
Beiter: 
~til4t 
se cong 
11283 

114349 


Cr ee 


11435 
wi lae7 
11383 
blows 
11739 
{i884 
11878 
etigee 
12826 
12189 
stom 
212248 
tee 
-f2o75 
12879 
12544 


cole 


ae 


SES 
Be frees 
9224 
Bye) 
9244 
7289 


7234 


IrgmM po YO ho Fo ro 


Pee er i eet 


2.9137 


mata 


DIA WE ANGLE 


8.287! 
9.2847 


vn 


ts, “Se e&% Sa wey “Sa ts Se ty SS em Ga cy, 
e s e es s e co e o e o e 
c 


‘Sa 
wee eee ee 


39979.9 TOY DEPTH: 29.2 
9.24 


WE DEPTH 
-289, 999 
linc vie 
=17hs oie 
192.998 
198,964 
“197.647 
~183.204 
~1899.049 
=l7euel) 
lite 
wiibae ae 
=) lees: 
169.144 
-146.924 
“Tae. /27 
=o25 400 
148,399 
159.247 
a Bebo et 
“134.819 
“131,724 
“149,899 
147.798 
145,747 
Ba ig 
“141.756 
139.798 
=1a76725 
alee ca 
oar ae 
=l51 3509 
a ae 
“128,985 
al coneae 
124, 288 
“128.456 
Si edva 
tle 
rd 
Sis weed 


—o5— 


1.9999 DE TOTAL ELONGATION 8.197 


HE REACH 

8.989 
{2.984 
24.178 
36.281 
48,394 
68.346 
72.648 
84.734 
Tonto 
199,971 


135.428 
145.688 
aa? 
169,992 
182.196 
196,486 
296.624 
218.848 
251,288 
243.318 
geveeo4 
267.815 
289.974 
Rage 

394,689 
316.886 
Serene 
341.459 
Sacer da 
366.956 
378.563 
398,678 
ad 
4toaod? 
427.649 
439.984 
452.528 
454,871 
477.928 
429, 008 
291.742 
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AHALYSTS PROGRAM OUTPUT 


SSN ATp 
11481 
11431 
1148! 
$1491 
11528 
elce? 
11589 
214489 
eiloot 
$1663 
ody 
pee) 
11736 
es 
L184 
11843 
stid72 
1199! 
Bees) 
Bee cus 
meee] 
212918 
212945 
12973 
mee 182 
RY 


41528. 417 
41929.912 
42967,496 
42387 .979 
42546. 758 
$2786.473 
43828.395 
43266. 213 
3386. 287 
43746. 28! 
43996. 434 
44226. 664 
44466.975 
44797.352 
44947, 389 
45188. 336 
#5428. 738 
42667,513 
49718. 300 
6151. 168 
46592. 83! 
4833. 954 
§5374, 929 
APG ie 
47526.204 
47597 .473 
47838. 746 
48689, 999 
42521.4998 
48562.96! 
43084433 
49948,374 
49297727 
49529. 4359 
39615.929 
oa254, 948 
29496, 844 
29738. 849 
29989.491 
ai222, 070 
31465, 169 
SMES ONL 
31949, 845 
s2191, 989 
52434, 349 
Sona. 76 
52919, 954 


mieoat St wwii 


SosOlara: 
wO484 557 
3647, 938 


556 


Ns 
£2996 
15864 


4 


e 
Poa C4 fo 


G4 roa mo oe 
O- ws = cy 


—J 


e :. ee & #4 
ee eee tee ee SE oe ee oY 


4 oa C4 ON C4 Oy Od ON CY OO 
oOo «yy Cr cn ON pe 

mu C4 O- OO fe GN 

tS qn ta ogo. = 


15884 
Pieuey 
14954 
$4185 
214393 
$4258 


s77 7 
0 $4335 


14497 
14482 
14537 
14632 
14797 
14782 
» 14857 


jet fat ta 


we 
no 
- 
x3 
a | 


“oO -O ~a 
“J “J oo 
cn -O cy ° 


hr ra IY or fr oP. 
s e e e 
4 
ond 
th 


oO oo OO Mm 


a Bit ft 


o 
o- 
ro 


ra 
« e * 
OH nwo WwW ow 
Oo wo -O 2. - 
an an 
nS 


937 


“oO ~3 
4 


J rah po PS oP PO on PS 
s s e e « 


ae it tm o 


toro 


e 
Oo co wg 
4) 
&S3 
oe | 


tate 


oo 3M ao © ow OI 


ed ee 2 OO oe Oe ee A ee Soe ee Oe Qe | 
e e « « s « « r 


8.43! 
3.1443 
8.1454 
§.1428 
3.1417 
3.1489 


78 


Re tsa 
« « e 
4% CA 
o~ ~~! 
cy = 


Eh ty “Se tse 
e«- ese  ® 


fd oy C4 Oa ON Oe OY Od 
om RY PI ON Be ote ON 


et tsa “Se esy 
Tm QO C4 -o CA po COO gn ro ~o O 


“Sr tS Sa esa 


tw rs EN GA 


hI 
co G- sy “So ~O -g Gm th 


tS. SSR tS Sy ty 
e & #« 


CF oe “I C4 OD ae 


Sa 
on 


an 
eS 


“Se Sy 
« « 


ee ce eo ee ee oe ee Oe Oe | 
C--4 be 
Oo on 


Sa tS 
e 


« 
bet pee OP 5 
OO 


er ne ee ee ee ee oe ee Oe ee oe ee 
~Q (Sa 32 


~O 
Oo ed Cre es OF igs CA wD 


“Sa eh “So th “Su es “Se te tha 
e so &® #« e 
~ 


Sa 
LJ 
-“~ oO co 


~~ 


Cr- 


« 
| cee oe coe aoe eee oe ed 


—s 


{So ey SSa cq “Sha tg “Sy 


Man <1 oO 
—- on -O C4 


= Niigeia 


“187.472 


-197,692 
Shean: 
~184.125 
-192,362 
169. 68 


4 CO rf 


a slg Sete 
Boy 
Bi oeeae 
SH ae) 
-72,994 
-79,486 
- 62,884 
SO gaat 
=tda) Oe 
“64,121 
“62.546 
68.779 
“39.416 
-37, 848 
“36,513 
“54.766 
Saree d 
eran el 
“39.171 
43.538 
-47,136 
45.626 
44,123 
“42.624 
sates! 
Seater 
aSuetts 
683 


i 
Cd ca & ES CH Cd 
Se re kd on 
« « 
~~. 
~be 
a. 


-86-— 


o14. 119 
o2h, 482 
308.069 
sul. 242 
265.652 
176.923 
388. 429 
689.823 
613.258 
625.642 
658.957 

39,401 
662.997 
Sicdaod 
od7 5474 
793,214 
712,859 


feap sid 


lide 


13/421 
Tho cky 
ope av 
824,862 
g57.3a1 
349.844 
$62.541 
SiesS 
387,343 
899.908 
Pierre 
924.998 
937.432 
eae 
Apeetot 
975. 893 
987.942 
1989. 384 
po12 5659 
1925.18! 
1937.733 
1956. 293 
1962.355 
1973.42 
$997,991 
1198565 
11i5.142 
Vio 4 
1138. 389 





ANALYSIS PROGRAM QUTFUT 


12468 
12495 


des 
pizcot 
esia 


12685 


oom 


Se awww 


euda?. 684 
34132.412 
24575. 285 
S4613.955 
Sete 
yui85.d44 
95946. 924 
Sade? G02 
Se rere 
36974.939 
30519, 287 
56562. 422 
26882. 689 
3784d.984 
ey ec 


ec? oo 


CA on OF on Cn 
fT oOo co oOo ds 


Sys 
14811 
16387 
018942 
17939 
slits 
17188 
17284 
aioe 
017445 
17499 
$7642 
hae 
STPRR’ 
2b 744 
10819 
18392 
18171 
018244 


r3G~ tS 


co co 0 co OO 


3 NSD 2 PD Pe PS 
* e 


mews td 


s 
co oO 


« 
ao 


eh po PO Po Ph PD PhO oP PQ 
e « . be e 
TO oo ta co cH aw (oO ow co a co 


Ca StS et oo ro An OO oC ON 


ans 
tym 


mn or 
om ron -O 


wng Sed my 


i | 
A © 
abts 


MN rh po ro 
eS e es 


thy “She esp fla egy CS 
ss» e&®© oe #&  #w «@ 


e 
| coe ene nena eee ce eee Oe ce Oe ed ed 


Ree by 6G 
aa | 


ee A Oy 


ro FO opt fed Cy Gd lp Re 


en eee coe a ne cee ol 
th Om t& 


cae on a cceel 
tS. ws, Gy 
PO cn -O 


~O 
oo 


The thy &Sy Oe Se ey, MSy ese Uy egy Ca 
- ##.e oo & #« se 8=6h la 


Sy be sy CSe ey 


‘si th 
s ro) 
— OF op  ~o -O 
CA op PO oon OO Po ON 


Sa ts 
« e 


me acne 
~27,922 


ais 7 inde 


=2Gnaid 
=e abo 
25.695 
Seed 
Sead 
Sl teoun 
coe 
~16.498 
we alee 
~15. 688 
~12. 296 
-19,899 
Suck 
Soa 
Snice 
SOedae 
So yy 
lene ths 
atee? 


-87- 


1139.878 
1143.499 


4eiue wi 


1244. 544 
1276.992 
1209 eau) 
1392.259 
1314.89! 
1ocieeoe 
{348.178 
Tog2ecee 
{365.474 
Poy saay 
1399.787 
1433,449 








ANALYSIS PROGRAN QUTPUT 


Smee RDIL281.0RT -- 


Pe elOsG DE TEN SE TAU GE DIA DE ANGLE 
85464 SI721.845 89852 2.976% §,2182 
20400 Sieiaioee cieee) 6 2.97b8 Bolo 
85098 52178.423 =. 99993, 2.9763) 9.2138 
ago! 32497.199 =. 188640029768) 9.2198 
89054 o2606.849 .i8i35 2.9706 9.298% 

NTE.) 92865.254 8.19297) 92.9753 8.2897 
83599 53894.652 18278 02.9758 = 9.2938 
Pub22 $5524.297 118549 82.9746 98.2815 
85644 Gogusnsos .tf42) 862.9745 68.1991 

o6h/ 33784.238 .19492 2.9749 43,1971 
95698 Suolangsce .itdo4 9) 2.7107 9.3998 
pear l2 $4244,986 1983500 2.97353500 9.1931 
BEY Ry) Gqeerctg ayer 269788 | 8.1912 
cise, parmo.soG | S778 2.7727 8 1895 
82773 Sivevemeom siveay. 2.0124 8. fais 
mecaes Borsceve)  bgeee = 2.9728 8, 1858 
85024 25498.898 8.19994 = 2.9717 = 8 1848 
22946 $2651.629 11866 2.9714 9.1824 
95969 Gdcoonvce 11138) 402,971) 8. 1887 
85891 Oheraenos) «6 11218002 279R BA L182 
span! Bhegieiao «=.11282 00294 
Pad Shemp. 276. . ide 2778) 9. 17h 
290957 So791.G74 = AL426 = 2.9898) = 9.1748 
Mey Segoe .itese 862.9695 838. 1732 
2699) Sjesaeded «tiatd 2.9682 9.1718 
86925 SASTyztl 9. ti6ss 892.9605 «=f. L798 
86845 Pees osOGmeligna 2.9885 8. Gy) 
6868 Bictetta «fae 06=—s 2 eel) 8 Lb? 
86989 oola7.976 .11869 2.9679 #8,1685 
sve. 1e S8421,064 .11932 2.9676 9.1632 
86132 38654.445 112884 02.9873 9, 1848 
#5153 oeda7.844 8.12877 892.9878 8. 1428 
95175 Demcice staiee «62.9057 8.1515 
86196 Biegewio. | .ie222 2.9663 8.1695 
86218 Giacusoro «i2dfee 2.96698 9.1593 
«95239 Etoeeesgo 8 .i2o0) «=. 2, PE oa LSOZ 
94241 46956.419 12848 2.9654 9.1371 
98292 49298.426 .12513 2.97651 9.1561 
Bb 594 49524.035 12585 2.9648 9.1959 
84523 AG7S8.745 12658 862.9643 8.1548 
88546 gov95.045 12751 2.9642 9.1538 
863567 41227,438 12884 802.9839 81521 
» 963288 SU4G1.9786 112875 2.9636 8.1511 
86415 42696.469 12949 82,9635 98, L582 
36431 AVoSte127 06 623822002. 76388 1492 
85432 42185.848 13992 2.7827 3.1483 
86473 42495.663 13168 2.9624 9.1474 


DE DEPTH 
-239. 888 
-197,434 
194.899 
=) 72ee2 
“180213 
197,46! 
“185.833 
aldeqocd 
~198, 26! 
“17 iega5 
ql) gaaee 
aoa 
“lia 
Seth [a3 
~166.495 
“164.278 
=1Geeprs 
ee a 
-137.736 
7150.074 
“193,449 
-131.565 
=| 492i 
-147, 299 
“145,144 
-143, 194 
=141.979 


=| i396 


~131.178 
129,249 
=litee ts 
-125,483 
125.996 
“121.628 
Senin ee 
117,985 
“116.935 
lai? 
Siler 
=li Ores 
199,735 
~184.956 
=i¥s, 172 
Slama ne 


—~88-— 


DE RESCH 
9, 389 
11.639 
23.288 
34,743 
46,411 
08. 289 
69.968 
81,698 
95.356 
185.841 
114,774 
128.494 
Novel 


is. ie] =U 
1o1.954 


137194 
178.954 
219.719 
222.498 
254.267 
246,958 
od) od 
Pehyorior 
281,43! 
EVes2ge 
580. B45 
316.969 
528.879 
543,084 
302, 304 
364, 148 
378, 297 
587,851 
579.679 
Ahad 
423.499 
447 137 
459.983 
479.982 
482,761 
494,544 
185.55! 
wL8. 422 
58.317 
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ANALYSIS PROGRAM QUTPUT 


@L59Q 


sa lw iw 


#5619 
86439 
» #8663 
86653 
1678) 
9572! 
95742 
86782 
» 86783 
95893 
86824 
pbd44 
78954 
95884 
» 86994 
96925 
88945 
86955 
#8985 
87995 
87925 
97945 
S700 
87994 
287194 
287124 
87144 
97165 
87133 
87293 
8/222 
87242 
87262 
$728! 


37438 
8745 

97474 
97493 


Guhug. dud 
42379, 527 
43195.979 
43349.787 
43975, 988 
45811,189 
44246,027 
44281.945 
44517,434 
44752.992 
44939,617 
49224. 393 
45469 .943 
45695.993 
45931. 786 
45167,719 
4£433.727 
#66379. 381 
45375. 934 
Wehlentcd 
47348379 
47534, 694 
47321 .947 
42057.467 
9293.94) 
48539.477 
43757963 
$7993,699 
49248395 
494771357 
49713.934 
49959,781 
J9197.878 


31513, 986 
31847, 329 
w2994.692 
eeoel.oa6 
Bocd?s 20) 
ayef publ 13) 
LSB a7 
Jo2/71, 689 
wos? 258 
Jo? 46. 924 
35944. 457 
Sieceelos 


24459, 892 


13241 
ASS14 
AB a) 
13469 
5 CHES 
15696 
13679 
leva 
» 13825 
13099 
lots 
14845 
14118 
1419! 
14285 


14330 


SBattraa 


14773 
1485! 
14925 
. 14993 
oe a 
15143 
eu. 
cde 
(asec 
16537 


idgis 


15586 
15669 
Ree 
£5897 
1508! 
S954 
16928 
14182 
016175 
124? 
wiGoeo 
16596 
16478 
16044 
»1ésta 
16692 
plevcc 
5 ito) 
{5913 


tr ves Pro RP ro roroormo Po op PS on PS op Por Po Sor Po Po og PO 
e * e bed e * e bed * e . rt) * s * 


. 
wo OO 2 
on 6. on 
os PD PD 
Oo IJ on 


ney) 


= 
O ~. 
q 
C. 


8.1466 
9.1457 
1449 
1449 
1432 
1424 
21417 
1439 
o.1491 
8.1594 
§.1397 
Sees 


Sa 


ee, “Sa thy “ha thy 
2 . 


Cad 
be | 
ho 


CN Cy C4 Cod ted CN Sed CA 


Sy oe FI pu CH Cy ke oon on of} 


Sa sy Sa chy 85s 8Sy 854 chy Sa egy ta 
s es es 2s s 
ao ho TS oO Po ON ep 6 On 


= 
C4 God 


a key 
tr- 


RG Sy esa %Sa ey Chr etsy 
s se . e 
“J sy OO .o 
CN ~O) «be thy 


t¥~ 


sa) PD sy PS omy Py PS 


ce i Oo a cee ee ee PO Lt a Pc 


3 4s | re Re cee ee Oe oe ee on oe De | oS oe 


tS. key 8S chy Sa tty Se oem 85s egy SS ogy 8S 
s es ® 2: » s zs. ese 8 
im oe ewe 


ro pa RO P 


Las] 
cn sa Gd oo Gd sy Ps 


oOo .0 


id 
tm ~yoms oop «OO ~o 
OM oy “sro OG es, 


—-89—-_ 


-191.636 
ieee 
=ooeae 
96.486 
94,692 
shuteo geile 
Sees 
-89,064 
-87.876 
86.196 
“84.525 
“82,9635 
“91.288 
SG e ys 
mito eae 
whi 28 
~74,579 
~73.954 
“71,446 
-59. 846 
69. 252 
66.666 
~65,387 
Sena 
=91 9748 
-§9. 339 
“38.907 
=o) eu! 
~55,752 


ww 8 


ra nls 
Seine 
=31, 178 
Bupa 
~48.146 
46.843 
45.146 
see otek 

42,168 
So ieeS 
“Siete 
Se I 
“34.829 
Aiticdiels 
BR 
Stati o 
oes 
727,691 
Rone ie 
~24,749 


vr 7San 
“23.028 


Soe 
346.924 
077.936 
3d? 859 
691.768 
615.699 
§25.4516 
637.545 
$49,478 
bé1.414 
Bl owed 
boseei7 
697,244 
799194 
721,147 
735,194 
740.300 
Poses 
fosve is 
788, 764 
792,949 
a94.716 
316.396 
928.983 
349,968 
852,956 
964.949 
376.944 
d82.343 
969.945 
912.958 
924,952 
956.396? 
748.084 
969.981 
BS 
934.944 
Tet od 
1388.99 
192!.939 
1833. 865 
1943, 162 
£957. 142 
199.185 
Meo ees! 
1993, 288 
Pivdesoc 
1117, 388 
1129, 443 
1141.583 
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ANALYSIS PROSRAN QUTPUT 


Pool S Gaao7moooe i070; geome 8.1163 21918 1143.63! 
Bed doz S4935.493 117961 892.9473 8 éh | | 28.595 AL ie 
87501 Beis Onecseiet tice 9) 207479 8. 1th7 17.898 Piet o9 
9735/8 Popliercomest cos 2.995) 8.itos -17.697 1291.943 
873597 go649.850 417282 2.9464 9.1149 -14. 399 ists 
87698 @ocey,g.6 .t/sc6 860 2,.9462 «9 T14G = 14,988 P220.776 
97827 26125.915 17439 2.9459 8.4142 -13.528 $2398.979 
27646 26365.803 17584 82.9406 8.2138 12.134 1208. 143 
87654 U6691.153 017978 82.7454 8.1154 0-18. 707 1262.208 
97683 Bedoveedemeisoge 2.74u1 9.188:  =9.922 1274.33? 
87792 37877.496 9.17725 02.9448 11270-8818 1286.45! 
97723 Siatg.cue, «t/088 2.9886 Oetl24 8 8=6-6.645 Pelened 
87748 Sucgseceae ot2d)a) 2.7485 8.1128 -37288 1319. 422 
27758 eiivfmee 17a Uke GIST 30821 2,22 
097777 ge058.485 18822 2.9438 g,4t13 -2.586 {334.824 
087798 Paceoweney «fosve 2.04385 8.1118 8 <te2ks 1545. 927 
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ANALYSIS PROGRAM OUTPUT 


WET DESIGN DATA DRY DESIGN DATS 
JR 233 
1143 9.24 $.1 
oH i 

$223 86322888 
8.999 
B.i83 869, 182 

VN ve Ri PH FT i WN VT RN FN FT X 
5.984 «14.679 =. 86399 8 8 § 3.245 14.644 = 92258 5 5 3 
Seeay (14e6S7)06=— SSeeC7989)0| 25.241) 12.884 = 3.288 «= 14.6530 091488 = 69S) 22.745 11.859 
5.914 14.694 = 84398 (esl) 0 2c.26) 24.178 See anita shoe) Go4ide 2257] 2) 2a. 2ee 
Beeermeeise 2) 864 7.858) 25.293 «636.2810 3.154 | 84.869 = 89388 Bed) 220079 eter as 
Qeevag 47970274 23517) (4865943899 14.876 = 888.1130 22.825 = 48. 41 
Geer 47id0| dicks 7,518 «625.341 8968.518 83.264 14.684 87iee 7.937 22.85 33.285 
Seoeoeia./2 de7ee 7.164 25.563 72.645 3.851 24.691 86339 7.764 22.3974 69.748 
Sgemete’co)6=6 7972907817) 25.096 «688786 | 29994697) b4de 74082 «= 22.897 = al. 68 
Byeereaerst 86778880 6.077) 023.487) 9894 28) 4.7840 8408507444 22,989 «= 93.558 
Beeeeetics)6=6 78288 )0=— 8741) «25.428 189.89 | 2737 M471 5688 7293022941 18S. 84) 
eee 6 fie | GGiL 625.848 t2l.255 0602987 A716 |= 288R 7146 22.982 «118.774 
2.743 14.747 78789 6.485 23,4867 133.429 2.878 14.721 81998 7.966 22.982 128.494 
2.797 14.751 76898 = 6.364 | 23.486 «145,688 2.02 14.727 g1198 6.87 23.991 149.221 
geaee 14,7560 752880247) 25.084 157.797 = 2,823 14.732 39-8 Gq. gus02 lolived 
Beef ol 6745980 GES4 23.522 169.992 «2.796 4737) |= 79889 =o B14 = 23.839 143.495 
Berometsres «6739880 6.8250 23.059) 192.198 2.77 34.742 78899 6.492 23.856 175.441 
fees. 69 = 73258 Dees conccore le, ane gee 14574) BIR Ob. S74 | 23,575 187,194 
eeg? 148.7730 72688 3.818 923.372 286.424 e/a | MARWoi 0 77ABG CG. 2h] 23.99 «(198.954 
eeveemette??)«=©6— 7898903719) 23.588 218.848 2.696 14.756 767968 6.151 23.386 219,719 
Bfecemmeryg: 4 7isee 865.624 «25.688 231.98 2.673 4.76 7688s 6.845 235.122 222.49 
feeeeeeean/0406=— 78786 06©=— 3.852 «(23.689 243.518 2.03 %4.764 73498 3.945 23.187 239.267 
2odi4 14,788 76198 3.442 23.034 255.564 2.428 14.759 74708 3.844 23,152 248.95 
Meemeece?))6=67G8S 0.556 |= 23.5449 «267.815 892.686 914.772 2 8=6A8G 3748 25,387 257.088 
eet? 14.795 889888 280.2720 «23.882 286.874 9 2.8d 4.778 ogee g0bo8 25.191) 247088 
Peeceeeewer?o 6 bU4e80 19) 25.676 292.538 2.564 14.779 72998 3.564 23.193 281.451 
2etos 14.081 67998 S.112 23.699 384.699 2.544 14.793 72588 3.477) 9 23.298 293.235 
Beeeesedor §6 AGB )|6OS 856 | 25.782 «316.886 2.524 14.786 71788 3.392 23,221 J85.545 
20079 14.887 60768 4.961 23.713 $29.17 2.584 14.789 71298 Sect 20.205 G1acge 
2.381 14.91 66338 894.889 23.728 344.459 2.485 14.793 97888 Gals 2oqc4o0 S2e0n09 
2.508 14,5815 65999 9.819 23.74 353.755 2.467 14.796 79108 3.153 23.258 340.594 
2.047 14.915 See yeid.  coupdemcooercs 2.999 14,799 69689 9.978 23.27 352.334 
2.531 $4,819 64999 4.685 23.764 378.363 2.431 14.992 49188 5.994 23.282 344,148 
eeesce «e990 A621 237s S98.b/6 0 6©=— 2.415 14.8850 B8bRE 4.95500 25,273 (375. 88? 
Meee teces §«6©6 GOR | 645E8 23.785 482.995 2.596 14.887 odie 4.864 23.584 397.85! 
2.205 14.925 63598 4.497 23.797) 413.919 Boe tee) | feet 46797) 23.515 377.677 
ge2od «(14.928 480983189) = 4.4380 25.888 427.849 2.363 14.813 7288 4.732 23.325 AL. S82 
Piegee itecs 6 627880 438k 25.819 «2459.984 «= 2547) 14.885 | fbbee 4.869 | 25.558 425.489 
Meese itscce «©6208 06=— 4524) «25.929 452.525 2.332 14.918 66208 4.687 23.348 435.271 
2.224 14.934 41998 4.27 25.859 464.671 2.3516 14.82 65899 4.547 23.556 447.337 

2.21 34.835 1508 4.216 25.95 477.8235 2.381 14.922 63589 4.488 23.565 459.808 
Smee l4dsco 6611824185 (25.859) «689.32 20287) 4825098844310 23.575 479.982 
2.183 14.84 $8788 4.414 23.9469 321.792 2.272 14.827 S4598 4.375 23.3584 482.741 
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BHALYSIS PROGRAM OUTPUT 


2.167 14.842 48389 4.965 23.879 St4.if 2.258 14,929 64198 4.521 23.593 494.644 
2.135 14.844 99999 4.917 23.888 926.492 2.244 {$4.831 63788 4.248 23.482 394.53! 
2.144 14.846 99489 Say 2oeuy) ) Jdenee eres aedoom | bused Angi? 23.41)° 318.422 
Ceereet-cnd 062880 S9248 623,986 Sol.242 0 2.217) «14.835 890962988 «= 4.166 225.42 938.517 
Bee ltsds)6©6CGemeS)0 3879) 95.915 «0985.63 902.284 57 25BE 4.117) «23.428 942.218 
2.197 $4,851 38889 3.935 23.924 376.923 2.191 14.839 62196 4.969 23.436 384.119 
Oeeeemtseovs «860220 S793 25.9350 88.4200 2.173) «= 14.84) 061888 = 4.823) 25,444 «366.826 
Peecewertecca «= 5/7998 203.7510 25.941 689.825 2.165 14.843 61498 3.977 23,492 377.936 
2.972 14.958 37099 Seyeeccstae oisees ) 2etdo) 14.045 GIdaNe 3.932 23.46 5 289.05 
2.685 14.859 37288 3.671 23.998 625.642 2.141 14.846 68788 3.869 25.468 681.74 
20849 814.859 36988 9.632 23.966 639.959 2.129 14.848 68408 3.846 23.474 613.69 
2.959 {$4,951 36688 3.994 23.974 658.481 2.118 14.85 68999 3.292 23.483 422.614 
Gece tface2 §86sbs89 = 3.556 025.982 662.987 2.195 14.851 39799 3.764 23.491 637.545 
2.817 14.3864 36899 Bhece conqgyordeece  2ogoa {4,d0o 9998 3.724. 25,490 649.478 
2.997 14.865 35788 3.485 23.997 687.774 2.884 $4,855 39489 3.435 23.595 661.414 
1,997 14.866 92439 AS 24.985 708.214 2.975 $4,956 38798 3.647 23.512 473.304 
fega, 0 $4,868 04090 5189346 |= 24.822 782.659 «= 262488 =: ag Ore ised ii) Gloue7) 
Meeeeco? 8 dgtcue 3.502 24.92 725.192 2,852 14.85% 99189 3.574 25.526 497.284 
mie 34.97 8634489 Soo coe eicfacoc erceUtcue taco. 8 Gold ~S.a8d 23.382 7996194 
Deomessc?2 0860880 3.S18 024.934 758.82 92.951 148620 | S789 LS 23.959 721.847 
1.949 14.373 54898 9.286 24.941 762.485 2.821 14.852 37283 3.469 23.545 733,184 
1.939 $4.974 93788 3.286 24.948 774.95 2.812 14.885 37889 3.438 23,552 743.985 
eoeeeAedss = su88 0 3.226 «624.855 787.421 2.982 14.966 36793 3.483 25.958 757.829 
1.972! 14.876 33288 3.196 24.862 799.997 1.992 14.867 36499 3.371 23.364 768.998 
1912 14.978 332799 3.167 24.869 812.577 1.983 14.968 56188 2.539 23.57 789.966 
1.994 14.9877 92799 3.139 24.975 824.9862 1.974 14.97 95999 3.588 23.577 792.94 
emcee teoo)06= 209803. ULL) «= o24.862 837,55! 1.965 $4.87) 33689 $.278 23.583 884.716 
1,887 $4.98) $2239 3.9984 24.999 849.944 1.954 14.872 85388 3.249 25.508 816.298 
1.878 84.982 32898 3.957 24.995 962.341 1.747) 14.875 9 35198 Sele § 20000 4uadeceas 
Meeeeeesees)«=— cee 03,95! 0 «(24.181 874.843 8 8091.938 $4,874 s488e 8 5.i9t 253.6 349.346 
Moe  ie.g84 06 S1Sde 03.885 | 24.187 987.548 1.929 i8.875 34689 3.163 23.686 852.83 
1.854 14.085 31588 2.98 24.114 999.958 {$.92! §4.876 34599 3.1356 23.611 464.249 
Meemeeerees |) «(911880029850 24.32 912.3572 1.913 14.878 34199 3.199 23.617 876.244 
Meccveeitece? 38996 2.931 24.125 924.89 1.994 14.879 33999 3.883 23.622 988.845 
83) 247288 «0s 8789020987) «24.152 937.412 1.896 14.88 35689 3.997 23.828 989,849 
1.924 14.989 28493 2,984 24.138 949.939 1.888 14.88! 35489 3.93! 23.6553 912.85 
1,816 $4.89 58298 2,04f 24.144 962.469 1.981 14.892 95299 3.896 25.839 924.3% 
1.889 4.891 58898 2.838 24.15 975.983 1.973 14.883 22998 2.992 23.644 936.08? 
1.992 14.891 49383 2.916 24.156 987.542 1.065 14.884 22798 2.9788 23.649 948.534 
1.795 $4,872 49698 2.774 24.161 1989,984 1.859 14.285 52588 2.934 23,654 969.991 
eccmetaec?s)§«©= 7488026775) 24.467 «1812.45 Gegeiaeeds sauce) 2,911 25.559 272572! 
Recmeeaesoe | 4929800 2.7520 24.173 1825.18! 1.843 {$4.996 92188 2.883 25,664 984.944 
774 14.995 89499800 2.7310 24.178 1937.735 891.886 14.887 31988 2.866 | 25.609 9996. 989 
1.767 14.896 48889 2.711 24.184 1958.293 1.828 14.8898 1783 2.9044 23.674 1998.798 
1.761 14.8796 48699 2.671 24.187 1962.855 1.921 14.989 $1588 2.822 23.478 1821.95 
Memeo’, «6«©6 8GG8Y 02,671) | 24,195 1873. 421 1.814 14,89 31398 2.891 23,403 1935.23 
1.748 14.999 48398 2.652 24.2 1907.99! 1.888 14,89) gig 2078 23,659 1945. 192 
1.741 14.879 48198 2.433 24.285 1189.565 %.88! 14.892 99999 2.757 25.693 19a/7.i4e 
1.732 14.8799 47993 2,614 24,211 1113.142 1.794 14.992 38789 2.739 23.697 1867, 185 
Py ss Pega 2rcse 9 24.216 $125,724 {1.768 14.893 $8588 2.717 25.782 1881.25! 
1.725 $4.781 47688 2.578 24.22) 1138.399 1.781 14.994  58o08 fae dd.to 1875.28 
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ANSLYSIS PROGRAM QUTPUT 


1.717 $4.98! 47493 2.06 24.226 1158.999 1.775 14.895  S81gg eee secee 1195. 552 
Peeie.col 0) «4/2990 2542 24.252 «(1165.49 861768 14.895 ©=—49908 = 2. 66E | 25.715 1117,386 
785) (14.983 «= 447888 2.525 0 24.257 1978.997 = t762 14.89%) = 9788264500 25,72 1129. 443 
1.699 {4.993 46998 2.588 24.242 11928.4687 1.756 14.3897 49688 2.424 23.724 1141.583 
Mecvomeeis.7o4 896788 2.492 24.247 $281.29 L753 14.998 49499 2.686 23.728 1153.56 
Weeceees.ge) 06 465990 2.475 0 24.252 (213.899 39st. 74t 14.89d = 492882 2.509 = 23.733 1165.63! 
1.4682 14.985 46499 2.459 24.257 1226.5! PV ccmelse 07) Mero eoue cdi te cdefos [177.699 
1.676 14.996 46298 2.443 24.261 1259.125 1.732 144.9 43996 2.554 23,741 1189.769 
1.8671 14.987 46195 2.427 24.266 251.743 1.726 14.9 48788 2.537 23.745 1281.843 
1,665 14.997 45999 2.412 24.27) 1264.3686 1.72 4.981 49588 POcome ecard) (fie. c19 
$66 °«(:14.992 8645888 2.597 | 24.276 1275.992 = 4 14.992 = R482 SBd 25.704 1225.998 
1,655 14.998 45698 2.382 24.281 1289.62! L789 14.962 48238 2.4889 23.758 1238.879 
1.649 14.999 45098 2.387 24.285 1382.255 1.783 14.993 48989 2.472 23.762 1258.13 
6446 14.910 40388 = 25520 24.27 1314.89! vot 0) Seen? 909 ecco 825.7 G0 00h lon. co 
Peover!) = 02990 2,538 924.295 $327,532 1.892 14.994 47788 2.441 23.77 1274.339 
boost 8614. 11 epee" 2.524 24.299 1548.175 1.687 14.985 47698 26425 23,774 1296.43! 
Megeeee 12.911 8 8=©644992 eee ectect a todzeced §s-502 914.905 eAvaag ate cde Vearoecce 
Meese si?) 478900 2296 24.588 1565. 4ie | LBP 986 |= 288 2.598 | 23.782 1318.622 
Memes tse 6 T4688 0 2.295) 24.585 1579.129 L.b7h 84.997 89497188 22.38) 23.783 1822.722 
b6t4 34.913 9 44598 02.2869) | 24.517 1399.787 = 8.866 14.9870 47888 2.387 0 25.789 1534.824 
el 14.9330 44398 02.256 |= 24.522 1493.449 «= bok 4.988 4bHBY 2.552 025.793 1546.929 
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